Astronomy & Astrophysics manuscript no. ms 


© ESO 2009 


October 8, 2009 





Evolution of the Early-Type Galaxy Fraction in Clusters since z = 

0.8* 

Luc Simard 1 , Douglas Clowe 2 , Vandana Desai 3 , Julianne J. Dalcanton 4 , Anja von der Linden 5 ' 6 , Bianca M. Poggianti 7 , 
Simon D. M. White 6 , Alfonso Aragon-Salamanca 8 , Gabriella De Lucia 6 9 , Claire Halliday 10 , Pascale Jablonka 11 , Bo 
Mil vang- Jensen 1213 , Roberto P. Saglia 14 , Roser Pello 15 , Gregory H. Rudnick 16 ' 17 , Dennis Zaritsky 18 

1 National Research Council of Canada, Herzberg Institute of Astrophysics, 5071 West Saanich Road, Victoria, British Columbia, 
Canada 

2 Ohio University, Department of Physics and Astronomy, Clippinger Lab 25 IB, Athens, OH 45701, USA 

3 California Institute of Technology, MS 320-47, Pasadena, CA 91 125, USA 

4 University of Washington, Department of Astronomy, Box 351580, Seattle, WA 98195-1580, USA 

5 Kavli Institute for Particle Astrophysics and Cosmology, PO Box 20450, MS 29, Stanford, CA 94309, USA 

6 Max-Planck-Institut fur Astrophysik, Karl-Schwarschild-Str., 1, Postfach 1317, D-85741 Garching 

7 INAF - Astronomical Observatory of Padova, Italy 

8 School of Physics and Astronomy, University of Nottingham, Nottingham, NG7 2RD, UK 

9 INAF - Astronomical Observatory of Trieste, via Tiepolo 11, 1-34143 Trieste, Italy 

10 INAF - Osservatorio Astronomico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy 

11 Observatoire de l'Universite de Geneve, Laboratoire d'Astrophysique de l'Ecole Polytechnique Federale de Lausanne (EPFL), 
1290 Sauverny, Switzerland 

12 Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark 

13 The Royal Library / Copenhagen University Library, Research Department, Box 2149, DK-1016 Copenhagen K, Denmark 

14 Max-Planck Institut fiir extraterrestrische Physik, Giessenbachstrasse D-85748 Garching, Germany 

15 Laboratoire d'Astrophysique de Toulouse-Tarbes, CNRS, Universite de Toulouse, 14 Avenue Edouard Belin, 31400-Toulouse, 
France 

16 The University of Kansas, Department of Physics and Astronomy, Malott room 1082, 1251 Wescoe Hall Drive, Lawrence, KS, 
66045, USA 

17 NOAO, 950 North Cherry Avenue, Tucson, AZ 85719, USA 

18 Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ, 85721 
Accepted for publication in A&A 

ABSTRACT 

We study the morphological content of a large sample of high-redshift clusters to determine its dependence on cluster mass and 
redshift. Quantitative morphologies are based on PSF-convolved, 2D bulge+disk decompositions of cluster and field galaxies on deep 
Very Large Telescope FORS2 images of eighteen, optically-selected galaxy clusters at 0.45 < z < 0.80 observed as part of the 
ESO Distant Cluster Survey ("EDisCS")- Morphological content is characterized by the early-type galaxy fraction f et , and early-type 
galaxies are objectively selected based on their bulge fraction and image smoothness. This quantitative selection is equivalent to 
selecting galaxies visually classified as E or SO. Changes in early-type fractions as a function of cluster velocity dispersion, redshift 
and star-formation activity are studied. A set of 158 clusters extracted from the Sloan Digital Sky Survey is analyzed exactly as the 
distant EDisCS sample to provide a robust local comparison. We also compare our results to a set of clusters from the Millennium 
Simulation. Our main results are: (1) The early-type fractions of the SDSS and EDisCS clusters exhibit no clear trend as a function of 
cluster velocity dispersion. (2) Mid-z EDisCS clusters around a = 500 km/s have f el ^ 0.5 whereas high-z EDisCS clusters have f et = 
0.4. This represents a ~25% increase over a time interval of 2 Gyrs. (3) There is a marked difference in the morphological content of 
EDisCS and SDSS clusters. None of the EDisCS clusters have early-type galaxy fractions greater than 0.6 whereas half of the SDSS 
clusters lie above this value. This difference is seen in clusters of all velocity dispersions. (4) There is a strong and clear correlation 
between morphology and star formation activity in SDSS and EDisCS clusters in the sense that decreasing fractions of [Oil] emitters 
are tracked by increasing early-type fractions. This correlation holds independent of cluster velocity dispersion and redshift even 
though the fraction of [Oil] emitters decreases from z ~ 0.8 to z ~ 0.06 in all environments. Our results pose an interesting challenge 
to structural transformation and star formation quenching processes that strongly depend on the global cluster environment (e.g., a 
dense ICM) and suggest that cluster membership may be of lesser importance than other variables in determining galaxy properties. 
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* Based on observations obtained in visitor and service modes at the 
ESO Very Large Telescope (VLT) as part of the Large Programme 
166.A-0162 (the ESO Distant Cluster Survey). Also based on ob- 
servations made with the NASA/ESA Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute, which is operated by 



the Association of Universities for Research in Astronomy, Inc., un- 
der NASA contract NAS 5-26555. These observations are associated 
with proposal 9476. Support for this proposal was provided by NASA 
through a grant from the Space Telescope Science Institute. 
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1. Introduction 

Our current paradigm for the origin of galax y morphologies rests 
upon hierarchical mass assembly (e.g., ISteinmetz & Navarrol 
2002), and many transformational processes are at work 
throughout the evolutionary histories of galaxies. Some de- 
termine the main structural traits (e.g., disk versus spheroid) 
while others only influence properties such as color and star- 
formation rates. Di sk galaxy collisions lead to the formation of 
elliptical galaxies ( Spitzer & Baade, 1951; Toomre & T oomrej, 
1971 iFarouki & Shaoirol 11982: Ne groponte & Whitel fl98H 



Barn es & Hernquistl Il992l 119961; iMihos & Hernquistl. [1996 
and the extreme example of this process is the build-up 
of the most massive galaxies in the Universe at the cores 
of galaxy clusters through the accretion of cluster mem- 
bers. Disks can also be transformed into spheroidals by tidal 
shocks as they are harass e d by t he cluster gr a vitational po - 
tential dFarouki & Shapirol Il98lt iMoore et all 1 1 9961 Il998l) . 
Harassment inflicts more damage to low luminosity galaxies be- 
cause of their slowly rising rotation curves and their low den- 
sity cores. Galaxies can be stripped of their internal gas and 
external suppl y through ram pressure exerted by the intraclus- 
ter m edium ( Gunn & G oto, 119721; iLarson. Tinslev. & Caldwelll 
1980; lOuilis. Moore. & Bowerl 120001) . and the result is a 
"quenching" (or "strangulation") of their star f ormation that 
leads to a rapid reddening of their colours (also see Martig ^eTaQ 
2009). The task of isolating observationally the effects of a given 
process has remained a major challenge to this day. 

Many processes affecting galaxy morphologies are clearly 
environmentally-driven, and galaxy clusters are therefore ideal 
laboratories in which to study all of them. The dynamical state of 
a cluster, which can be observationally characterized by measur- 
ing mass and substructures, should be related to its morpholog- 
ical content. For example, the number of interactions/collisions 
suffered by a given galaxy should depend on local number den- 
sity and the time it has spent within the cluster. Dynamically 
young clusters with a high degree of subclustering should con- 
tain large numbers of galaxies that are infalling for the first time. 
More massive clusters will contain more galaxies, but they will 
also have high er galaxy-galaxy r elative velocities that may im- 
pede merging dLubin et all 120021) . Spheroidal/elliptical galaxies 
will preferentially be formed in environments where the bal- 
ance between number density and velocity dispersions is op- 
timal, but it is still not clear where this optimal balance lies. 



ity dispersion dRood et al.Lll972llDresslei 


1984: Carlberg et al. 


1997; iTranetal.. 19991 Borganietal. 


19991 Lubinetal. 


2002), through weak-lensing shear (Kaiser & Squires, 1993 



Schneider & Seitz, 1995; Hoekstra, Franx, & Kuiiken, 2000; 



Clowe et all 120061) or th rough analys is of their hot X-ray emit- 



ting atmospheres (e.g., lAllenl 1 1998b . and it will be used here 
as the main independent variable against which morphological 
content will be studied. 

The morphological content of high-redshift clusters is most 
often characteri zed by the fraction /e+so of early-type galax- 



Fasano et al., 2000; van Dokkum et al., 2001; Lubin et al., 2002; 


Holden et all l2004t ISmifh et all 


2005: iPostman et all 120051 


Desai et al., 2007; Poggianti et all 


2009b). The bulk of the data 



available so far is based on visual classification. "Early-type" 
galaxies are defined in terms of visual classifications as galax- 
ies with E or SO Hubble types. A compilation of ear l y-type 
fractions taken from the literature (Ivan Dokkum et all [2000) 
shows a dramatic increase of the early-type fractions as a 



function of decreasing redshift from values around 0.4-0.5 at 
Z ~ 1 to values around 0.8 in the local Universe. However, 
the inter pretation of this trend is not ent i rely clear as oth- 
ers (e . g.. iDressler et all 119971 ; iFasano et all 120001 ; iDesai et aTl 
120071; iPoggianti et a l.. 2009bJ) have reported that the fraction 
of E's remains unchanged as a function of redshift and that 
the observed changes in early-type fractions are entirely due 
to the SO cluster populations. SO populations w ere observed 



to grow at the expense of the spiral populati on (Smith et al. 



2005; IPostman et all 12005; Mora n et all 120071; IP oggianti et al- 
2009bJ although others (e.g.. iHolden et all 120091) have argued 



for no e volution in the relati ve fr action of ellipticals an d SOs with 
redshift. ISmith et all d2005l) and lPostman et alj d2005l) show that 
the evolution of /e+so is in fact a function of both lookback time 
(redshift) and projected galaxy density. They find /e+so stays 
constant at 0.4 over the range 1 < ti 00 kback < 8 Gyr for projected 
galaxy densities S < 10 Mpc~ 2 . For high density environments 
(I = 1000 Mpc~ 2 ), f E+s0 decreases from 0.9 to 0.7. At fixed 
lookback time, /e+so varies by a factor of 1.8 from low to high 
densities at ti 00 kback — 8 Gyr and by a factor of 2.3 at ti 00 kback = 1 
Gyr. The difference between low and high density environments 
thus increases with decreasing lookback time. Both studies in- 
dicate that the transition between low and high densities occurs 
at 0.67?2oo (^2oo is the projected radius delimiting a sphere with 
interior mean density 200 ti mes the critical dens ity at the clus- 
ter redshift, see Equation [TJ. Postman et all (120051) also find that 
/e+so does not change with cluster velocity dispersion for mas- 
sive clusters (<x > 800 km/s). The data for one of their clusters 
also suggest that /e+so decreases for lower mass systems. This 
trend would be consistent with observations of /e+so in groups 
that show a strong trend of dec reasing /e+so versus decreasing <x 
(Zabludoff & Mulchaey, 1998). Finally, f E+S o seems to correlate 
with cluster X-ray luminosity at the 2-3cr level dPostman et all 
l2005h . 

R ecent works on stel l ar mass-selected cluster galaxy sam- 
ples (IHolden et all 120071; Ivan der Wei et all 120071) paint a dif- 
ferent picture. The fractions of E+SO galaxies in clusters, 
groups and the field do not appear to have changed signifi- 
cantly from z ~ 0.8 to z ~ 0.03 for galaxies with masses 
greater than 4xlO IO M . The mass-selected early-type fraction 
remains around 90% in dense environments (E > 500 gal 
Mpc~ 2 ) and 45% in groups and the field. These results show 
that the morphology-density relation of galaxies more massive 
than 0.5M* has changed little since z ~ 0.8 and that the trend 
in morphological evolution seen in luminosity-selected samples 
must be due to lower mass gala xies. This is in agree ment with 
iDe Lucia et aU d2004l 120071) and lRudnick etaTI (120091) who have 
shown the importance of lower mass (i.e., fainter) galaxies to the 
evolution of the color-magnitude relation and of the luminos- 
ity function versus redshift. Another interesting result has come 
from attempts to disentangle age, morphology and environment 
in the Abell 901/902 supercluster d Wolf et al I 120071: lLane etall 
2007). Local environment appears to be more important to 
galaxy morphology than global cluster properties, and while 
the expected morphology-density and age-morphology relations 
have been observed, there is no evidence for a morphology- 
density relation at a fixed age. The time since infall within the 
cluster environment and not density might thus be the more fun- 
damental parameter dictating the morphology of cluster galax- 
ies. 

A number of efforts have been made on the t heoretical side 
to mo del the morphological content of clusters. Diaferi o~et all 
(2001) used a model in which the morphologies of cluster galax- 
ies are solely determined by their merger histories. A merger be- 
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tween two similar mass galaxies produces a bulge, and a new 
disk may form through the subsequent cooling of gas. Bulge- 
dominated galaxies are in fact formed by mergers in smaller 
groups that are later accreted by clusters. Based on their model, 
they reach the following conclusions: (1) the fraction of bulge- 
dominated galaxies inside the virial radius should depend on the 
mass of the cluster, and it should show a pronounced peak for 
clusters with mass of 3 x 10 14 M Q followed by a decline for 
larger cluster masses. (2) The fraction of bulge-dominated galax- 
ies should be independent of redshift for clusters of fixed mass, 
and (3) the dependence of morphology on cluste r mass should 
be str onger at high redshift than at low redshift. lLanzoni et al.l 
(2005) use the GALICS semi-analytical models and find that 
early-type fractions strongly depend on galaxy luminosity rather 
than cluster mass. By selecting a brighter subsample of galax- 
ies from their simulations, they find a higher fraction of ellip- 
ticals irrespective of the cluster mass in which these galaxies 
reside. This trend is particularly noticeable in their high-density 
environments. Observations and these earlier models clearly do 
not agree in important areas, and a comparison between them 
would clearly benefit from a larger cluster sample size. More re- 
cently, the Millennium Simulation (MS; ISpringel etail [2005) 
has provided the highest resolution model thus far of a large 
(0.125 Gpc 3 ), representative volume of the Universe. Improved 
tracking o f dark matter structure and new semi-analytical pre- 
scriptions (iDe Lucia & Blaizoi, 120071) allow the evolution of the 
galaxy population to be followed with higher fi delity and bet- 
ter sta tistics than in the otherwise similar work of Di aferio et al.l 
(120011) . We will use cluster catalogues from the MS later in this 
paper for comparison with our observational data. 

Our understanding of high-redshift cluster galaxy popu- 
lations in terms of their evolution as a function of redshift 
and their cluster-to-cluster variations has been hampered by 
the lack of comprehensive multi-wavelength (optical, near- 
infrared and X-ray) imaging and spectroscopic studies of 
large, homogeneously-selected samples o f clusters. Many efforts 
are underway to improve sample sizes ([G onzalez et al.. 2001; 
iGladders & Yeel 120051: IWillis et al.Ll2005llPostman et alll2005l) . 
One of these efforts is the Europea n Southern Obser vatory 
Distant Cluster Survey ("EDisCS"; IWhite et al.L 120051) . The 
EDisCS survey is an ESO large programme aimed at the study 
of a sample of eighteen optically-selected clusters over the red- 
shift range 0.5-0.8. It makes use of the FORS2 spectrograph on 
the Very Large Telescope for optical imaging and spectroscopy 
and of the SOFI imaging spectrograph on the New Technology 
Telescope (NTT) for near-inf rared imaging. A number of p apers 



on star formation in clusters (Poggianti et al., 2006, 2009a) and 



the assembly of the cluster red sequence dDe Lucia et al.L 120041 
l2007tlSanchez-Blazquez et1dll2009tlRudnick et al.Ll2009l) have 
been so far published from these data. In addition to the core 
VLT/NTT observations, a wealth of ancillary data are also be- 
ing collected. A 80-orbit program for the Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope was devoted to 
the /-band imaging of our ten highest-redshift clusters. Details 
of the HST /ACS observation s and visual galaxy classifications 
are given in lDesai et alj d2007l) and the frequenc y and properties 
of galaxy bars is studied in lBarazza et al.l (12009). X-ray observa- 
tions with the XMM-Ne wton satellite of thr ee EDisCS clusters 
have been published in lJohnson et alj ([2006) with more clusters 
being obser ved. H-alpha obser vations of three clusters have been 
published in lFinn et alj {2005) with more clusters also being ob- 
served. Finally, the analysis of Spitzer/IRAC observations of all 
EDisCS clusters is in progress (Finn et al., in preparation). 



This paper presents the early-type galaxy fractions of 
EDisCS clusters as a function of cluster velocity dispersion, 
redshift and star-formation activity. A set of local clusters ex- 
tracted from the Sloan Digital Sky Survey (SDSS) is used 
as a comparison sample. Early-type fractions were measured 
from two-dimensional bulge+disk decompositions on deep, op- 
tical VLT/FORS2 and HST/ACS images of spectroscopically- 
confirmed cluster member galaxies. Section [2] describes the 
EDisCS cluster sample selection and the imaging data. SectionfS] 
describes the procedure used to perform bulge+disk decomposi- 
tions on SDSS, VLT/FORS2 and HST/ACS images. Section H 
presents early-type fractions for the EDisCS clusters with a de- 
tailed comparison between visual and quantitative morphologies 
and between HST- and VLT-derived early-type fractions. It also 
includes early-type fractions for the SDSS clusters. Changes 
in EDisCS early-type fractions as a function of cluster veloc- 
ity dispersion, redshift and star-formation activity are studied 
in Section [5] Finally, Sections [6] and [7] discuss our results and 
their implications for the morphological content of clusters. The 
set of cosmological parameters used throughout this paper is 
(H ,n m ,n A ) = (70,0.3,0.7). 



2. Data 

2.1. Sample Selection and VLT/FORS2 Optical Imaging 

The sample selection and optical/near-infrared i maging data for 
the EDisCS survey are de scribed in details in iGonzalez et al.l 
d2002l) . IWhite et al] d2005l) (optical photometry) and Aragon- 
Salamanca et al. (near-IR photometry; in preparation). 
Ph otometric r e dshif ts for the EDisCS clusters are presented 
in iPello et all (120091) . and cluster velocity dispersions mea- 
sured from weak-lensing mass reconstructions are given in 
IClowe et alj d2006[). Spectrosc opy for the EDisCS clusters is de - 
tailed in lHallidav et all d2004) and lMilvang- Jensen etai] (l2008). 
Clusters in the EDisCS sample were drawn from the Las 
Campanas Distant Clu ster Survey (LCDCS) candidate catalog 
dGonzalez et all 1200 ll) . Candidate selection was constrained by 
published LCDCS redshift and surface brightness estimates. 
Candidates were selected to be among the highest surface bright- 
ness detections at each redshift in an attempt to recover some 
of the most massive clusters at each epoch. Using the esti- 
mated contamination rate for the LCDCS of ~ 30%, we tar- 
geted thirty candidates in the redshift range 0.5-0.8 for snapshot 
VLT/FORS2 imaging in an effort to obtain twenty (10 at z ~ 0.5 
and 10 at z ~ 0.8) confirmed clusters. 

The z ~ 0.5 candidates were observed for 20 minutes in 
each of Ib and Vb, and the z ~ 0.8 candidates were observed 
for 20 minutes in each of 7 B and R s? . These filters are the stan- 
dard FORS2 o nes. Vb and 7 B are close approximations to the 
iBesselll (Il990l) photometric system while the /? sp is a special fil- 
ter for FORS2. Final cluster candidates for deeper VLT imag- 
ing were selected o n the basis of colo r and surface density of 
galaxies on the sky (IWhite et all 120051) . The image quality on 
th e final stacked ima ges ranged from 074 to 078. As described 
in IWhite et aU (120051) . deep spectroscopy was not obtained for 
two cluster candidates (1122.9-1136 and 1238.5-1144), and we 
therefore did not include them here. The main characteris- 
tics (positions, redshifts, velocity dispersions and radii) of the 
EDisCS cluster sample used in this paper are given in Table Q] 
7^200 is the projected radius delimiting a sphere with interior 
mean density 200 times the critical density at the cluster redshift, 
and it is used throughout this paper as an important fiducial ra- 
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dius. R200 values in Table Q] were calculated using the equation 



^200 - 1-73 



cr 1 
lOOOkm/s -y/Q A + Q m (l + Z ) 



(1) 



where /i 100 = Ho / 100 and cr c i uster is the cluster velocity disper- 
sion measured using spectroscopically-confirm ed cluster mem- 
bers (ICarlberg et all 11997b iFinn et all l2005h . Cluster masses 
were calculated using the equation: 



M c , = 1.2 x 10 



i— , 

\1000km/s/ Ji 



1 M Q (2) 



as in lFinn et alj d2005l) . 

In practice, the redshift distributions of high-z and the mid-z 
samples partly overlap as can be seen from Table Q] 

2.2. VLT Spectroscopy and Cluster Membership 

We use only spectroscopically-confirmed cluster members to 
calculate our cluster early-type fractions. Deep multislit spec- 
troscopy of the EDisCS was obtained with the FORS2 spectro- 
graph on VLT. Spectra of > 100 galaxies per cluster field were 
obtained with typical exposure times of two and four hours for 
the mid-z and high-z samples respectively. Spectroscopic tar- 
gets were selected from /-band catalogues. This corresponds 
to rest-frame ~ 5000 + 400 A at the redshifts of the EDisCS 
clusters. Conservative rejection criteria based on photometric 
redshifts were used in the selection of spectroscopic targets to 
reject a significant fraction of non-members while retaining a 
spectroscopic sample of cluster galaxies equivalent to a purely 
/-band selected one. We verified a posteriori t hat these crite- 
ria excluded at most 1 % of the clu ster galaxies (Hal lidav et al.L 
2004; Milvang ~Jensen et aTl l2008). The spectroscopic selection, 
observations a n d spe ctroscop ic catalogs are presented in detail in 
lHallidav et all d2004l) and iMilvang- Jensen et all (120081) . As de- 
scribed in Hal lidav et al.1 d2004l) . cluster redshifts and velocity 
dispersions were iteratively calculated using a biweight scale es- 
timator for robustness. Cluster members were defined as galaxies 
with redshifts within the range Zduster ± 3cr c / !wfer where Zduster is 
the median redshift of all cluster members. 

2.3. HST/ACS Imaging 

In addition to our ground-based imaging, a 80-orbit program 
(GO 9476, PI: Dalcanton) for the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST) was devoted to the 
/-band imaging of our ten highest- redshift cluster fi elds. Details 
of these observations are given in lDesai et al.l(l2007l) . Briefly, the 
HST observations were designed to coincide as closely as pos- 
sible with the coverage of the ground-based optical imaging and 
spectroscopy, within guide star constraints. The VLT/FORS2 
images cover a 6'5 x 6'5 region around each cluster, with the 
cluster center displaced by l'from the center of the region. For 
reference, the ACS WFC has a field of view of roughly 3'.5 x 3 '5. 
Balancing scientific motives for going deep over the entire spec- 
troscopic field against a limited number of available orbits, we 
tiled each 6'5 x 6'.5 field in four 1 -orbit pointings overlapping 
one additional deep 4-orbit pointing on the cluster center. The 
resulting exposure time per pixel was 2040 seconds except for 
the central 3'.5 x 3'.5, which had an exposure time per pixel of 



10200 seconds. The deep central pointing probes to lower sur- 
face brightness, fainter magnitudes, and larger galactic radii in 
the region of the cluster containing the most galaxies. All expo- 
sures were taken under LOW SKY conditions to maximize our 
surface brightness sensitivity. An image mosaic was created for 
each cluster using the CALACS/Multidrizzle pipeline, and the 
final sampling of the multidrizzled image mosaics was 07045. 
This is the "native" ACS image sampling, and it was chosen 
to avoid potential aliasing problems that might have been intro- 
duced by a finer multidrizzle sampling given our limited dither 
pattern in the cluster outskirts. Clusters with HST imaging are 
identified by a "h" in Table Q] 

3. Quantitative Galaxy Morphology 

3. 1 . Source Detection and Extraction 

The source catalogs and segmentation images for the EDisCS 
clusters were created using the SExtractor (" Source Extractor") 
galax y photometry package version 2.2.2 dBertin & Arnoutsl 
1996). The SExtractor source detection was run on the combined 
deep FORS2 images in "two-image" mode using the I-band im- 
age as the reference detection image for all the other passbands. 
The detection threshold was 1 .5cr^, and the required minimum 
object area above that threshold was 4 pixels. The convolution 
kernel was a 7x7 Gaussian kernel with a FWHM of 3.0 pix- 
els. No star/galaxy separation based on the SExtractor "stellar- 
ity" index was attempted. Every source was fit with a bulge+disk 
model, and unresolved sources such as stars could easily be iden- 
tified as output models with zero half-light radius. 

As SExtractor performs source detection and photometry, it 
is able to deblend sources using flux multi-thresholding. This 
deblending technique works well in the presence of saddle 
points in the light profiles between objects. Each SExtractor pre- 
deblending "object" consists of all the pixels above the detection 
threshold that are spatially connected to one another. This group 
of pixels may or may not include several real objects. The multi- 
thresholding algorithm assigns the pixels between two adjacent 
objects and below the separation threshold based on a probabil- 
ity calculated from bivariate Gaussian fits to the two objects. No 
assumption is made regarding the shape of the objects in this sta- 
tistical deblending technique. We used a value for the SExtractor 
deblending parameter DEBLEND MINCONT of 0.0005. This 
value is subjective, and it was found through visual inspection 
of several EDisCS cluster images to provide good object sep- 
aration. Even though the value of DEBLEND MINCONT was 
determined subjectively, it provides an unequivocal definition of 
an object in the EDisCS catalogs. It was only determined once, 
and the same value of DEBLEND MINCONT was consistently 
used for all EDisCS cluster images as well as for all the reliabil- 
ity tests of Section [J.2.5l 

3.2. Two-Dimensional Bulge+Disk Decompositions 

This work uses GIM2D (Galaxy I Mage 2D) ver s ion 3 .2, a 
2D decomposition fitting program (Sim ard et all l2002f) . to 
measure the structural parameters of galaxies on the EDisCS 
VLT/FORS2 and HST/ACS images. GIM2D is an IRAFQ'SPP 
package written to perform detailed bulge+disk surface bright- 
ness profile decompositions of low signal-to-noise (S/N) images 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. Main characteristics of the EDisCS cluster sample: IDs, positions, redshifts, number of spectroscopically-confirmed mem- 
bers, velocity dispersions and radii. Clusters with HST imaging are identified by the superscript "h" in their ID. 



Mid-z clusters 



ID 


RA a 


DEC 


u 

Z h 


Age of Universe 


/V c 

1 v mem 


,1 

cr d 




* AT f 




(2000.0) 


(2000.0) 




(xt ) 




(km/s) 


(Mpc) 


(1O 15 M ) 


(1) 


(a) 


(3) 


(4) 


(5) 


(0) 


(H\ 
(l) 


(8) 


(9) 


1018.8-1211 


10:18:46.8 


-12:11:53 


0.4716 


0.654 


33 


474 + " 
-37 


0.91 


0.142 


1059.1-1253 


10:59:07.1 


-12:53:15 


0.4550 


0.663 


41 


517 + 71 
165 +fj 


1.00 


0.186 


1119.3-1130 


11:19:16.7 


-11:30:29 


0.5491 


0.615 


21 


0.30 
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11 Cluster BCG Coordinates (12000) 
b Cluster redshift measured from EDisCS spectroscopy 
c Number of cluster members confirmed by EDisCS spectroscopy 
d Cluster velocity dispersion measured from EDisCS spectroscopy 
c From equation [Tj 
f From equation|2] 



of distant galaxies in a fully automated way. GIM2D is publicly 
available, and it has been used extensively in a wide range of 
different projects so far. 

3.2.1. Fitting Model 

The fitting model used for the two-dimensional bulge+disk de- 
compositions of ED isCS galaxies is the same as the one used by 
Sima rd et al.l (|2002). It consists of a "bulge" component with a 
de Vaucouleurs profile and of an exponential "disk" component. 
We put "bulge" and "disk" between quotes to emphasize that this 
conventional nomenclature does does not say anything about the 
internal kinematics of the components. The presence of a "disk" 
component does not necessarily imply the presence of an actual 
disk because many dynamically hot systems also have simple 
exponential profiles. The fitting model had ten free parameters: 
the total galaxy flux F, the bulge fraction B/T (= for pure disk 
systems), the bulge semi-major axis effective radius r e , the bulge 
ellipticity e (e = 1 - b/a, b = semi-minor axis, a = semi-major 
axis), the bulge position angle of the major axis (pi, on the image 
(clockwise, y-axis = 0), the disk semi-major axis exponential 
scale length rj (also denoted h in the literature), the disk inclina- 
tion i (face-on = 0), the disk position angle </></ on the image, the 
subpixel dx and dy offsets of the model center with respect to 
the input science image center. The sky background is not a free 
parameter of the fits (see Section [3.2.3l ). The Sersic index for the 



bulge profile is fixed at a value of n — 4 (i.e., the de Vaucouleurs 
profile value). The position angles 0j and were not forced to 
be equal for two reasons: (1) a large difference between these 
position angles is a signature of strongly barred galaxies, and (2) 
some observed galaxies do have bona fide bulges that are not 
quite aligned with the disk position angle. 

The smooth bulge+disk model used here is obviously a sim- 
ple approximation. After all, many real galaxies will exhibit 
more than two structural components such as nuclear sources, 
bars, spiral arms and HII regions. Even in the presence of only 
a bulge and a disk, the ellipticity and/or the position angles of 
these components might be functions of galactocentric distance. 
The bulge+disk model is a trade-off between a reasonable num- 
ber of fitting parameters and a meaningful decomposition of dis- 
tant galaxy images. No non-parametric or parametric quantita- 
tive classification system is perfect. Any classification system 
will suffer from biases inherent to its basic definition. However, 
provided a given quantitative system is clearly defined before its 
use, its results will be readily reproducible in their successes and 
failure by other investigators. 

The exact shape of bulge profiles remains under debate (e.g., 
Balc ells et al.L [2003, and references therein). Locally, there is ev- 
idence that the bulges of late-type spiral galaxies may be better 
fit by an n — 1 profile, whereas bright ellipticals and t he bulges 
of ea rly-type spiral g alaxies follow an n = 4 profile dde Jond , 
ll996UCourteau et aU[T99llAndredakislll998l) . Local late-type 
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galaxies with n — 1 bulges have B/T < 0.1 dde JongL [1996). 
Since such bulges contain only 10% of the total galaxy light, low 
signal-to-noise measurements of late-type high-redshift galaxies 
make it very difficult, if not impossible, to determine the Sersic 
index of distant bulges even with the spatial resolution of the 
Hubble Space Telescope as demonstrated by an extensive set 
of tests on HST im ages of the high-redshift cluster CL1358+62 
dTran et all l2003h . On the other hand, n is more important for 
bulge-dominated galaxies, and n = 4 is the expected value based 
on local early-type galaxies. Knowing that bright ellipticals and 
the bulges of early-type spirals are well-fit by a de Vaucouleurs 
profile, a n — 4 bulge profile was therefore adopted as the canon- 
ical bulge fitting model here for the sake of continuity across the 
full range of morphological types. 

3.2.2. Fitting Regions 

GIM2D disk+bulge decompositions are performed on thumb- 
nail (or "postage stamp") images extracted around the objects 
detected by SExtractor rather than on the entire science image 
itself. The area of the thumbnail images is given by the isophotal 
area of the object. Here, all thumbnails were chosen to have an 
area 5 times larger than the 1.5cr Mg isophotal area. Each thumb- 
nail is a square image with sides of length y5 x isophotaljarea. 
The first thumbnail is extracted from the science image itself, 
and the local background calculated by SExtractor is subtracted 
from it so that it should have a background mean level close to 
zero. The second thumbnail is extracted from the SExtractor seg- 
mentation image. The GIM2D decompositions were performed 
on all pixels flagged as object or background in the SExtractor 
segmentation image. Object areas in the segmentation image are 
sharply delineated by the location of the isophote corresponding 
to the detection threshold because SExtractor considers all pixels 
below this threshold to be background pixels. However, precious 
information on the outer parts of the galaxy profile may be con- 
tained in the pixels below that threshold, and fits should there- 
fore not be restricted only to object pixels to avoid throwing that 
information away. Pixels belonging to objects in the neighbor- 
hood of the primary object being fit are masked out of the fitting 
area using the SExtractor segmentation image. The flux from 
the primary object that would have been in those masked areas 
in the absence of neighbors is nonetheless properly included in 
the magnitude measurements given in this paper because mag- 
nitudes were obtained by integrating the best-fit models over all 
pixels. 

3.2.3. Sky Background Level Measurements 

Special care must be paid to the determination of the local sky 
background level b and dispersion <Xbk g as sky errors are the 
dominant source of systematic errors in bulge+disk decompo- 
sitions of distant galaxies. As an example, overestimating the 
background sky level will lead to underestimates of the galaxy 
total flux, half-light radius and bulge fraction as a result of strong 
parameter covariances. Even though the SExtractor local back- 
ground was subtracted from each galaxy thumbnail image, an 
additional (residual) background estimate db was computed and 
used by GIM2D to correct for any systematic error in the initial 
SExtractor sky level estimate. In order to compute db, GIM2D 
used all the pixels in the science thumbnail image flagged as 
background pixels (flag value of zero) in the SExtractor seg- 
mentation image. GIM2D further pruned this sample of back- 
ground pixels by excluding any background pixel that is closer 



than five pixels (l'.'O for the pixel sampling of the FORS2 de- 
tectors) from any (primary or neighboring) object pixels. This 
buffer zone ensures that the flux from all SExtracted objects in 
the image below all the l.Scbkg isophotes does not significantly 
bias the mean background level upwards and artificially inflate 
o~bkg- A minimum of 7500 sky pixels was imposed on the area 
of the sky region. In cases where the number of sky pixels in 
the input science thumbnail image was insufficient, the origi- 
nal science image was searched for the 7500 sky pixels nearest 
to the object. For the EDisCS fits, background parameters were 
re-calculated with GIM2D before fitting, and the residual back- 
ground levels db were then frozen to their recalculated values for 
the bulge+disk fits. 

3.2.4. Point-Spread-Functions 

The shape of the Point-Spread-Function (PSF) on the 
VLT/FORS2 and HST/ACS images varies significantly as a 
function of position, and these variations must be taken into ac- 
count when point-spread-functions for the bulge+disk decompo- 
sitions are generated. For both sets of images, we used the stand- 
alone ve r sion o f the stellar photometry program DAOPHOT II 
(Stetson, 1987) to construct spatially-varying PSF models for 
the EDisCS cluster images. For each cluster and for each pass- 
band, we selected "clean", point sources (detection flag of zero 
and stellarity index of 0.8 or greater) from the SExtractor source 
catalog. The positions of these point sources were fed to the 
DAOPHOT routine PSF to be modelled as the sum of a Gaussian 
core and of an empirical look-up table representing corrections 
from the best-fitting Gaussian to the actual observed values. 
Both the gaussian core parameters and the look-up table were 
allowed to vary linearly as a function of x and y positions on the 
image. Finally, the PSF model was used to create a PSF at the 
position of each galaxy to be fit. The PSF images were 2'.'5 on a 
side to provide good dynamical range for the fits. 

3.2.5. Reliability Tests 

Following the same procedure as in [Simard et al. (2002]), we 
performed an extensive set of simulations to test the reliabil- 
ity of our sky background estimates and of the best-fit param- 
eter values recovered through bulge+disk fits on both sets of im- 
ages. 2000 smooth galaxy image models were created with struc- 
tural parameters uniformly generated at random in the following 
ranges: 20.0 < I < 25.0, 0.0 < B/T < 1.0, < r e < 1070, 
0.0 < e < 0.7, < r d < 1070, and 0° < i < 85°. The bulge 
Sersic index was held fixed at n — 4 for all models. Both bulge 
and disk position angles were fixed to 90° for all simulations, 
and the bulge and disk sizes were uniformly generated in the log 
of the size ranges above. Each simulation was convolved with a 
PSF computed from one of the images with a FWHM typical of 
the VLT/FORS2 (~ a/8) and HST/ACS (~ 0705) observations. 
The same PSF was used in both creating and analyzing the sim- 
ulations, so the results will not include any error in the structural 
parameters due to PSF mismatch. Poisson deviates were used to 
add photon noise due to galaxy flux into the simulations. The 
noisy images were then embedded in a 20" x 20"section of one 
of the real /-band images to provide a real background for the 
simulations. In addition to sky photon noise and detector read- 
out noise, the real background noise includes brightness fluctu- 
ations of very faint galaxies below the detection threshold. This 
procedure thus yields realistic errors that include the effect of 
sky errors. The simulations were SExtracted exactly in the same 
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way as real EDisCS sources (see Section [3~TT i. Science and seg- 
mentation thumbnails extracted from the simulations were ana- 
lyzed with GIM2D following exactly the same steps as for the 
real galaxies (see Section [X2| >. 

Figures Q] and [2] show maps of errors on the galaxy total 
magnitude /, galaxy intrinsic half-light radius r^i and galaxy 
bulge fraction B/T for the VLT/FORS2 images. The left-hand 
panels show the mean parameter errors as a function of input 
galaxy magnitude and size, and the right-hand panels show the 
lcr parameter random error as a function of input galaxy mag- 
nitude and size. The lower number in each cell is the num- 
ber of simulated galaxies created for that cell. Most system- 
atic errors are directly related to surface brightness as magni- 
tudes and sizes of low surface brightness sources are inherently 
harder to measure. This fact is borne out by the trends in the er- 
rors shown in Figure Q] Decreasing surface brightness follows a 
line going from the lower left-hand corners to the upper right- 
hand ones. The top panels of Figures Q] show that systematic 
errors on / start to become significant (A/ 0.2) fainter than 
/ = 22.5. Systematic errors on log r/,/ also increases signifi- 
cantly beyond this magnitude. It is important to note that / = 
22.5 is significantly fainter by about 2 mag than the galaxies 
that will be used to compute cluster early-type galaxy fractions 
in Section 14.31 so these galaxy fractions should be unaffected. 
Figure [2] shows that systematic errors on B/T are smallest over 
the region / < 22.5, -0.5 < log ru < 0.3 where most of the 
real EDisCS galaxies actually lie. As mentioned above, our re- 
liability tests do not include the effects of PSF mismatch errors 
because we used the same PSF for creating simulated images 
and for their analysis. However, we were able to check that these 
errors were not significant because we fitted both galaxies and 
stars on our real VLT/FORS2 images. The measured intrinsic 
radii of the stars clustered at zero, and this would not have been 
the case should PSF mismatch errors have been important. 

4. Early-Type Galaxy Fractions 

4. 1 . Definition and Comparison with Galaxy Visual 
Classifications 

The bulk of the previous work on the morphological content 
of high-redshift clusters is based on the visual classification 
of galaxies, and this section compares visual and quantita- 
tive morphological classification. Visual classifications for 9200 
ga laxies in EDisCS clusters with HST images are presented 
in IDesai et al. J 2007I) . As shown by previous works (|lm et all 
2002t iMcIntosh et all l200l iTran et all I2003L iBlakeslee et all 
2006), quantitative and visual morphologies can be best linked 
together by focussing on three structural parameters: bulge frac- 
tion B/T, image smoothness 5 and bulge ellipticity e. The image 
smoothness, 5 , is defined as: 

S=R t +Ra (3) 

where Rj and Ra are defined in Equation 1 1 of Simard et all 
(2002). These two indices quantify the amount of light in sym- 
metric and asymmetric residuals from the fitting model respec- 
tively, and they are expressed as a fraction of the total galaxy 
model flux. 5 is typically measured inside a radius that is a 
multiple of the galaxy half-light radius. Using our HST/ACS 
measurements, we found no differences between image smooth- 
ness within one and two galaxy half-light radii. We therefore 
use image smoothness inside two half-light radii (and denote 
it 52 hereafter) because it is more reliably measured on the 
VLT/FORS2 images with their lower spatial resolution. We can 



choose selection criteria on B/T, 5 and e that yield the best 
match to the visual classifications, and the particular choices are 
not important as long as the same selection criteria are applied 
to both local and high-redshift clusters. 

We divide the visually-classified EDisCS into T = — 5 (E), 
-2 (SO), 1 (SO/a) and "others" (T > 1). Using our HST/ACS 
structural parameter measurements, we find that E and SO galax- 
ies have similar B/T distribution with the SO distribution being 
skewed towards slightly lower B/T, but e distributions are dif- 
ferent. It is therefore possible to differentiate between E and SO 
galaxies on the basis of these two parameters. SO and SO/a galax- 
ies have similar e distributions but different B/T and 5 distribu- 
tions. Given that the bulge ellipticity e cannot be reliably mea- 
sured on the VLT/FORS2 images, we restrict on selection crite- 
ria to B/T and 52. Figure[3]shows 52 versus B/T for the four vi- 
sual types of galaxies. 52 can take on small negative values du e 
to statistical background subtraction terms (Sim ard et all l2002). 
The optimal choice of limits on B/T and 5 2 for our definition of 
early-type fraction is driven by the need to maximize the num- 
ber of E/S0 galaxies selected while minimizing the contamina- 
tion from Sa-Irr galaxies. After several iterations, we settled on 
B/T > 0.35 and 52 < 0.075 as our definition of an early-type 
galaxy. These limits are very similar to those used in previou s 
studies dim et al. l l2002tlMcIntosh et alll200llTran et alll2003l) . 
With these criteria, our quantitative selection can be translated 
into visual classification terms as 



U =(0.69Afe + 0.71/V S0 + 0.35jV SO / fl + 0.04N Sa - Ir r)/N total (4) 

The coefficients in Equation [4] give th e completeness of the 
quantitative classification in terms of the lDesai et alj d2007l) vi- 
sual classes. For example, the adopted B/T and 52 cuts would 
select 69% of the galaxies visually classified by IDesai et ail 
(120071) as E's, 71% of their SO's and so on. As mentioned 
earlier, E's and SO's cannot be distinguished using only B/T 
and 52. Equation [ 4] is t o be compared to the prescription of 
Ivan Dokkum et allfcoOOh : 



1 

Set =(Ne + Ne/so + Nso + 



T S0/a)/N, 



total 



(5) 



where N total is the number of galaxies with My < -20. 

It is impossible to recover all the galaxies visually classified 
as early-types because a visual early-type does not necessarily 
imply a r 1 ^ 4 profile. Indeed, many early-type gal axies such as 
dwarf ellipticals have si mple exponential profiles dLin & F aber. 
1983; Kormendy, 1985), and we have verified through isophote 
tracing that many galaxies visually classified as early-types and 
missed by our selection criteria do have radial surface bright- 
ness profiles that are exponential and thus consistent with their 
measured low B/T values. 

Given N tota i galaxies brighter than an absolute magnitude 
limit Myjhn inside a clustercentric radius R max of which N et are 
early-types galaxies, we actually calculate the early-type galaxy 
fraction by finding the median of the binomial probability distri- 
bution 



p{x)dx - 



AW 



-x N «(l-xf 



(6) 



N et \{N tot ai-N et )\ 

and we integrate Equation [6] to calculate the lower and upper 
bounds of the corresponding 68% confidence interval. In the 
limit of large Ntotai an d N e , (not always true for the current clus- 
ter sample), this converges to the same symmetric error bars as 
would be obtained from the propagation of gaussian errors. 
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Fig. 1. Two-dimensional maps of GIM2D systematic and random galaxy magnitude and half-light radius errors from 2000 
VLT/FORS2 image simulations. Top left-hand panel: Systematic error on recovered galaxy total magnitude I rec as a function of 
input galaxy log half-light radius r n ij nput in arcseconds and input galaxy total magnitude Ii„ pul . The top number in each cell is the 
mean magnitude error (I rec - I input), and the bottom number is the number of simulations created in that cell. Top right-hand panel: 
1 cr random error on I rec (o~(I rec — ^m/jw/)) a. function of log r fiij n p U t 

and Impui- Bottom left-hand panel: Systematic error on recovered 
galaxy intrinsic log half-light radius r/,/ rec as a function of input galaxy log half-light radius r n ij np ut in arcseconds and input galaxy 
total magnitude Ii„ put . The top number in each cell is the mean log radius error (log rhi, rec - log r n i,input), and the bottom number is 
the number of simulations created in that cell. Top right-hand panel: lcr random error on log r n [ rec (cr(log r/,; rec - log r„i input )) as a 
function of log r UMput and I input . 



4.2. HST-Based Fractions 

For each EDisCS cluster with HST/ACS imaging, we have com- 
puted the fraction of early-type galaxies using our quantitative 
HST/ACS morphologies (B/T > 0.35 and and 52 < 0.075). We 
used only spectroscopically-confirmed members brighter than 
an absolute V-band magnitude Myjim- We varied %,„ as a func- 
tion of redshift from -20.5 at z = 0.8 to -20.1 at z = 0.4 to ac- 
count for passive evolution. This choice of My.u m was made to 
be fully consistent with previous work dPoggianti et all [2006) 
although it may not be strictly the best choice for late-type 
galaxy populations. Our results did not appear to be sensitive to 
variations in Myn m at the l evel o f a few tens of a magnitude. 
Following Poggian ti et alj (120061) . our early-type galaxy frac- 
tions were also computed by weighting each galaxy according 
to the incompleteness of the spectroscopic catalog. This incom- 
pleteness depends on both galaxy magnitude and clustercentric 



position. Incompleteness as a function of magnitude was com- 
puted by dividing the number of galaxies in the spectroscopic 
catalog in a given magnitude bin by the number of galaxies in 
the parent photometric catalog in the same bin. We used 0.5 mag 
bins here. Incompleteness due to the geometrical effects comes 
from the finite number of slitlets per sky area, and the increasing 
surface density of galaxies on the sky closer to the cluster cen- 
ters. Geometric incompleteness is field dependent as it depends 
on cluster richness, and we thus computed this incompleteness 
on a field-by-field basis. We also used four radial bins out to 7?2oo 
with a bin width of 0.25^2oo- 

The raw and incompleteness-corrected HST-based early- 
type galaxy fractions are given in Table [2] for a maximum clus- 
tercentric radius R et of 0.6/?200 (columns 4 and 5) and /?200 
(columns 9 and 10). Most of the corrected fractions do not sig- 
nificantly differ from the raw ones because our spectroscopic 
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Fig. 2. Two-dimensional maps of GIM2D systematic and random galaxy bulge fraction errors from 2000 VLT/FORS2 image simu- 
lations. Top left-hand panel: Systematic error on recovered galaxy bulge fraction (B/T) rec as a function of input galaxy log half-light 
radius rhi,i np ut in arcseconds and input galaxy total magnitude Ij nput . The top number in each cell is the mean bulge fraction error 
((B/T) rec - (B/T)inp Ut ), and the bottom number is the number of simulations created in that cell. Top right-hand panel: Icr random 
error on (B/T) rec (cr((B/T) rec - (B/T) input )) as a function of log r MMput and I input . 



sample is essentially complete down to / < 23 {My ~ -20 
at z = 0.8), and we used multiple masks on dense clusters 
to improve the spatial sampling of our spectroscopic sample. 
As a comparison, Table [2] also gives early -type galaxy frac- 
tions measured from visual classifications bv lDesai et alj (12007) 
(Columns 6 and 7). They should be compared with values in 
Column 5 because cluster galaxy samples selected using pho- 
tometric redshifts are de facto free from the magnitude and ge- 
ometric incompleteness of our spectroscopic sample. Another 
important caveat is that they were computed using two different 
ways to isolate cluster members (photometric redshift and statis- 
tical background subtraction), and they are thus not restricted to 
spectroscopically-confirmed members. Nonetheless, the agree- 
ment between fractions measured from visual and quantitative 
classifications is remarkably good. The largest disagreement is 
for 1 138.2-1 133, but even this case can be considered marginal 
as it is not quite 2cr. 



4.3. VLT- versus HST-based Fractions 

Quantitative morphologies measured from HST images are 
more robust than those measured fro m ground-based images 
(Section |3~2~5l and lSimard et al.ld2002h ). Figure|4]shows a direct 
galaxy-by-galaxy comparison between bulge fraction and im- 
age smoothness measurements from HST/ ACS and VLT/FORS2 
images. This comparison includes spectroscopically-confirmed 
member galaxies from all clusters with HST imaging that are 
brighter than Myjim and within a clustercentric radius of O.6.R200 
to take into account the effect of crowding. For a given galaxy, 
the agreement between the two sets of measurements will ob- 
viously depend on its apparent luminosity and size. The over- 
all agreement is reasonably good. The scatt er in the bulge frac- 
tion plot is consistent with ctb/tacs ~ 0.1 ( Si mard et all [2002) 
and ctb/t,vlt ~ 0.25 (Figure|2]) added in quadrature, but the fact 
that completely independent segmentation images were used for 
the HST and VLT morphological measurements also contributes 
significantly to this scatter. Indeed, this scatter would be smaller 
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Table 2. Early-Type Galaxy Fractions Based on HST/ACS Imaging 



ID R et < 0.6^200 Ret S ^200 





^*clus 


Net 


fet,raw 


fet,corr 


fEISO.phz 


JElSQMg 


N ctus 


Ne, 


fet,raw 


fei,corr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1216.8-1201 


45 


23 


0.51 ±0.07 


0.55+0.07 


0.47+0.06 


0.54+0.06 


57 


25 


0.44+0.06 


0.42+0.06 


1040.7-1156 


9 


4 


0.45±0.15 


0.45+0.15 


0.53+0.17 


0.37+0.17 


13 


5 


0.40+0.12 


0.33+0.12 


1054.4-1146 


18 


9 


0.50±0.11 


0.29+0.10 


0.28+0.09 


0.24+0.09 


26 


11 


0.43+0.10 


0.28+0.09 


1054.7-1245 


11 


8 


0.70±0.13 


0.46+0.14 


0.44+0.16 


0.57+0.13 


19 


10 


0.52+0.11 


0.38+0.11 


1232.5-1250 


48 


28 


0.58±0.07 


0.58+0.07 


0.60+0.06 


0.45+0.06 


51 


28 


0.55+0.07 


0.47+0.07 


1037.9-1243 


8 


1 


0.18±0.12 


0.18+0.12 


0.13+0.09 


0.08+0.07 


8 


1 


0.18+0.12 


0.18+0.12 


1103.7-1245b 


2 





0.21±0.20 


0.21+0.20 


0.47+0.20 


0.21+0.20 


4 





0.13+0.13 


0.13+0.13 


1354.2-1231 


8 


5 


0.61+0.15 


0.39+0.15 


0.35+0.14 


0.44+0.19 


12 


5 


0.42+0.13 


0.28+0.11 


1138.2-1133 


22 


4 


0.20±0.08 


0.20+0.08 


0.37+0.09 


0.50+0.14 


24 


6 


0.26+0.08 


0.34+0.09 



a From Table 14 of Desai et al. ( 2007) 
b From Table 16 of besai et all (2007) 
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Fig. 3. Image smoothness parameter S2 versus bulge fraction 
B/T for different visual types. The galaxies selected by our 
quantitative early-type galaxy criteria {B/T > 0.35 and S2 < 
0.075) are enclosed in the area delimited by dashed lines. 



if only uncrowded galaxies (as indicated by the SExtractor pho- 
tometry flag) on the VLT images had been plotted here. For the 
image smoothness plot, there is a correlation between SlfORsi 
and S 2 acs > but it is not one-to-one. S 2acs values increase faster 
than S2fous2- This is expected as PSF blurring will be more sig- 
nificant on the ground-based images, and 52 measurements are 
not corrected for PSF effects. Part of the scatter is again due to 
the use of independent segmentation images. 

The inclusion of clusters with only VLT/FORS2 imaging al- 
lows us to extend our analysis to nine additional clusters - an 
important consideration given that we seek to probe cluster-to- 
cluster variations in morphological content. We therefore need to 
show that we measure consistent early-type fractions for clusters 
with overlapping ACS and FORS2 images. The problem boils 
down to finding the set of limits on B/TfORsi and S2fORsi that 
yield FORS2 early-type fractions in agreement with the ACS 
fractions obtained with B/Tacs ^ 0.35 and S2acs ^ 0.075 
when the same galaxies are used for both FORS2 and ACS. 
For each cluster, we used all spectroscopically-confirmed cluster 
members brighter than Myjim an d within a clustercentric radius 
of /?200- No corrections for incompleteness were applied here 
as these corrections would be identical for both cases. We went 
through many manual iterations until we found satisfactory lim- 
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Fig. 5. Comparison between early-type galaxy fractions for clus- 
ters with overlapping VLT and HST imaging. VLT/FORS2 
and HST/ACS early-type galaxy fractions were computed us- 
ing galaxies with B/Tforsi ^ 0.40 and S2fORsi < 0.05 and 
B/T acs > 0.35 and S2 AC s < 0.075 respectively. The ACS and 
FORS2 f el values plotted here are listed in column 4 of Table [2] 
and column 4 of Table [3] Dashed line is the one-to-one line. 



its on B/Tforsi and S2forsi- We found FORS2 fractions to be 
in very good agreement with the ACS ones for B/Tforsi > 0.40 
and S2fors2 ^ 0.05 (Figure 0. This agreement is especially 
good if one considers the fact that we performed our FORS2 
and ACS bulge+disk decompositions completely independently 
from one another, i.e., we did not attempt to use the same 
SExtractor segmentation map for both FORS2 and HST images. 
The limit on B/Tfors2 is slightly higher than the one on B/T acs 
because lower spatial resolution typically leads to a small over- 
estimate of the bulge fraction. Similarly, the limit on S2fors2 
needs to be more stringent than on S2acs to select the same 
galaxies as they will look smoother on the FORS2 images due 
to lower resolution. 
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Fig. 4. Direct galaxy-by-galaxy comparison between bulge fraction (left-hand panel) and image smoothness (right-hand panel) 
measurements from HST/ACS and VLT/FORS2 images. Filled circles are galaxies classified as early-type on both ACS and VLT 
images, asterisks are galaxies classified as early-type only on the VLT images, pluses are galaxies classified as early-type only on 
the ACS images, and open circles are galaxies not classified as early-type on either ACS or VLT images, The dashed lines show the 
cuts used for the definition of an early-type galaxy as discussed in Sections [4. 1 I and l4~3l 



Following the procedure described in Section l4~2l we com- 
puted early-type galaxy fraction for all eighteen clusters using 
galaxies on our FORS2 images with B/Tforsi ^ 0.40 and 
S2fORsi ^ 0.05. The results are shown in Table [3] The same 
incompleteness corrections as in Section l4~2l were applied here 
as well. The errors on the early-type galaxy fractions in the table 
do not include errors on /?200 due to correlated errors on cluster 
cr. We hereafter use our VLT/FORS2 early-type fractions for all 
EDisCS clusters for the sake of uniformity. 

4.4. Local Clusters 

The Sloan Digital Sky Survey (SDSS: lAbazaiian et al.Ll2009h of- 
fers by far the best, "local" (z < 0.1) baseline for a comparison 
of early-type galaxy fractions between local and high-redshift 
clusters. Clusters similar in mass to EDisCS clusters can be se- 
lected from spectroscopic SDSS data, and galaxy morphologies 
can be measured using GIM2D from SDSS images. We there- 
fore used SDSS-selected clusters here to construct a local base- 
line as nearly free of systematics as currently possible given the 
available data. 

We use t h e sam ple of SDSS clusters defined in 
Ivon der Linden! (120071) . T he basis of this c luster sample is 
the C4 cluster catalogue dMiller et al.L 120051) . and we briefly 
recapitulate here how the von der Linden et al. sample was 
selected. Their primary aim was to find the galaxy closest to the 
deepest point of the potential well of a cluster. In order to insure 
that the clusters would span a large angular extent compared to 
the minimum distance of 55 arcsec between fibers, the sample 
was restricted to redshifts z < 0.1. This first cut resulted in an 
initial sample of 833 clusters. A combination of clustercentric 
distance, galaxy concentration and colour cuts was used to 



identify brightest cluster galaxies (BCGs) for these clusters. For 
cases where the same BCG was identified for more than one 
cluster, only the cluster with the density peak was retained, and 
the others were deemed to be substructures. This cut rejected 
101 clusters. Refined velocity dispersion and virial radii were 
then computed through an iterative process of velocity cuts. This 
process failed for 55 clusters, and these were also rejected. All 
remaining clusters were then visually inspected. An additional 
set of 35 clusters were rejected at this point as being in the 
infall regions of other clusters, and another 17 clusters were 
discarded because they had less than three galaxies within 3<x 
of the cluster redshift and l/?200 of its center. T his brought the 
total o f SDSS clusters down to 625. Following iPoggianti et al] 
(12006b . we applied a final redshift cut to keep clusters in the 
range 0.04 < z < 0.085. The lower limit reduces fiber aperture 
effects, and the upper limit minimizes incompleteness in galaxy 
absolute magnitude. Our final SDSS comparison sample thus 
has 439 clusters. 

Given that we are interested in probing galaxy properties 
as a function of environment, it is important to ensure that the 
SDSS and EDisCS samples both cover the same range of envi- 
ronments. We therefore selected a subsample of SDSS clusters 
with a velocity dispersion distribution matching the EDisCS dis- 
tribution. This match was done by adding SDSS clusters to the 
subsample one at a time and keeping only those that maintained 
the EDisCS-SDSS two-sample Kolmogorov-Smirnov probabil- 
ity above 50%. This is the probability of the maximum difference 
between the normalized cumulative distributions of the EDisCS 
and SDSS samples. It means that even if the two sampls were 
selected at random from the same underlying distribution, they 
would differ by more than the two observed samples more than 
half the time. This probability threshold thus yields a SDSS sub- 
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Table 3. Early-Type Galaxy Fractions Based on VLT/FORS2 Imaging 



ID 






R e , < 0.6R 


200 




Ret < R 


200 






N il 


fet,raw 


fet,corr 


N il 


fet,raw 


fet,corr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


1018.8-1211 


-20.2 


18 


0.60+0.11 


0.55+0.11 


20 


0.59+0.11 


0.55+0.11 


1059.1-1253 


-20.2 


20 


0.64±0.10 


0.59±0.10 


28 


0.64±0.09 


0.60±0.09 


1119.3-1130 


-20.3 


6 


0.64±0.17 


0.77±0.15 


9 


0.64±0.14 


0.74±0.15 


1202.7-1224 


-20.1 


11 


0.54±0.14 


0.54±0.14 


13 


0.60±0.13 


0.67±0.13 


1232.5-1250 


-20.2 


48 


0.48+0.07 


0.60±0.07 


51 


0.49+0.07 


0.62±0.07 


1301.7-1139 


-20.2 


17 


0.53±0.11 


0.53±0.11 


28 


0.43±0.09 


0.40±0.09 


1353.0-1137 


-20.3 


9 


0.55±0.15 


0.55±0.15 


17 


0.42±0.11 


0.36+0.11 


1411.1-1148 


-20.2 


15 


0.47±0.12 


0.53±0.12 


16 


0.44±0.12 


0.44±0.11 


1420.3-1236 


-20.2 


4 


0.50±0.20 


0.50+0.20 


7 


0.56±0.17 


0.44±0.17 


1037.9-1243 


-20.3 


8 


0.18±0.12 


0.18±0.12 


8 


0.18±0.12 


0.18+0.12 


1040.7-1156 


-20.4 


9 


0.36+0.14 


0.36±0.14 


13 


0.46+0.13 


0.46±0.13 


1054.4-1146 


-20.4 


18 


0.40+0.11 


0.34±0.11 


26 


0.39±0.09 


0.35±0.09 


1054.7-1245 


-20.5 


11 


0.62±0.13 


0.38±0.13 


19 


0.52±0.11 


0.43±0.11 


1103.7-1245b 


-20.4 


2 


0.21+0.20 


0.21+0.20 


3 


0.39+0.22 


0.39±0.22 


1138.2-1133 


-20.2 


22 


0.50±0.10 


0.46+0.10 


24 


0.54±0.10 


0.54±0.10 


1216.8-1201 


-20.5 


45 


0.47±0.07 


0.53±0.07 


57 


0.46+0.06 


0.46±0.06 


1227.9-1138 


-20.3 


9 


0.26±0.13 


0.16+0.11 


11 


0.30+0.12 


0.22±0.11 


1354.2-1231 


-20.5 


8 


0.50+0.16 


0.39+0.15 


12 


0.35+0.12 


0.28+0.11 



a Number of cluster members brighter than M vljm inside R, 



sample that is very well-matched to the EDisCS clusters. The re- 
sulting subsample (referred to as "SDSS-C4" hereafter) includes 
158 clusters, and these clusters are listed in Table|4] 

We ran GIM2D on SDSS Data Release Seven (DR7; 
Abazajian et al. , 2009) u-, g-, r- and /-band images of objects 
in the magnitude range 14 < r petros i anzCorr < 17.77 with a 
galaxy spectrum (i.e., with field SpecClass = 2 in database ta- 
ble SpecPhoto). Bulge+disk decompositions were successfully 
obtained for 674,693 galaxies (Simard, in preparation). GIM2D 
morphologies for galaxies in our matched SDSS-C4 clusters 
were extracted from this large morphological database to com- 
pute early-type fractions. There are two sources of incomplete- 
ness that must be taken into account here. The first one is in- 
completeness versus magnitude. We denote this spectroscopic 
completeness function as C mag (m) here, and we compute it 
around each cluster position by taking the ratio of the number 
of galaxies in the spectroscopic SDSS catalog (database table 
SpecPhoto) to the number of galaxies in the photometric SDSS 
catalog (database table PhotoPrimary) as a function of Petrosian 
r magnitude. Galaxies around a given position on the sky were 
extracted from the database using the SDSS "fGetNearbyOb- 
jEq" function. The second source of incompleteness comes from 
the spatial sampling of the SDSS fibers on the sky. Fibers can- 
not be placed closer than 55" from one another. This means that 
regions with a higher surface density of targets could not be sam- 
pled as completely as regions in the global field. The net result 
for SDSS clusters is a decrease in spectroscopic sampling as a 
function of decreasing clustercentric distance R. We can map the 
spectroscopic completeness versus R by computing the ratio of 
galaxies in the spectroscopic and photometric SDSS catalogs as 
a function of R. We denote this geometrical completeness func- 
tion as C geom (R) here. Ideally, C geom (R) should be computed for 
each cluster because it will depend on cluster richness and ap- 
parent size (and thus indirectly on redshift). However, in prac- 
tice, there are not enough galaxies in a single cluster to yield 
C ge0 m(R) with acceptable error bars. So, we opted for averaging 
clusters with the same redshifts and velocity dispersions to com- 
pute C geom (R). We divided the cluster list of Table [4]into three 
cluster groups: (1) z < 0.06, (2) z > 0.06, cr < 800 km/s, and 



(3) z > 0.06, cr > 800 km/s. The weight W, pec (m, R) in the spec- 
troscopic catalog of a galaxy with a r'-band magnitude m at a 
clustercentric R is thus given by the product r —. r 1 , and 

the completeness-weighted early-fraction of a SDSS cluster is 
then simply: 

2j W spec (m h Rj) 
ie[Af v <-19.8, 
B/7>0.35, 




Wspec(mi,Rd 



i6[Af v <-19.8, 
«<««] 

In terms of spatial resolution, the ACS, SDSS and FORS2 
images have sampling of 0.68 kpc/FWHM at z = 0.8 (a/09 
FWHM in i), 1.87 kpc/FWHM at z = 0.07 (l'/4 FWHM in 
g) and 4.5 kpc/FWHM at z = 0.8 (a/6 FWHM in /) respec- 
tively. Even though the sampling of the ACS and FORS2 im- 
ages differs by a factor of seven, their limits on B/T and S2 for 
the computation of consistent early-type galaxy fractions were 
quite similar. This is an indication of the robustness of our mea- 
sured structural parameters over this range of spatial resolutions. 
For the sake of simplicity, we therefore adopt the ACS limits 
((B/T) SD ss, g > 0.35 and S2 SDS s,g < 0.075) for our SDSS early- 
type galaxy fractions rather than use yet another set of limits. 
We can further test these limits on the catalogue of visually 
classified galaxies f rom the SDSS North Equatorial Region of 
Fukugita et al. (2007). This catalogue contains Hubble T-type vi- 
sual classifications for 2253 galaxies down to a magnitude limit 
of r — 16. If we apply our limits on (B/T)sDss,g an d S2$Dss,g 
to galaxies in this catalogue, then we find that the coefficients 
of the SDSS-to-visual equivalent of Equation [4] would be 0.88, 
0.68, 0. 14, and 0.014 respectively. Early-type SDSS galaxies are 
therefore quantitatively selected with an "efficiency" compara- 
ble to our selection from the ACS images. 

The raw fractions of [Oil] emitters for the 158 SDSS-C4 
clusters were calculated by directly querying the SDSS database 
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Fig. 6. Comparison between fractions of [ QUI emitters com- 
puted using emission-line measurements from lBrinchmann et al.l 
(2004) and the DR7 release. Filled and open circles are clusters 
with cr > 600 km/s and cr < 600 km/s respectively. 



table SpecLine for the [Oil] 3727 and [Oil] 3730 equivalent 
widths for each confirmed cluster member, adding them together 
and correcting them to rest-frame by dividing by (1+z). The 
corrected [Oil] fractions were then computed following exactly 
the same calculations (and using the same weights, the same 
luminosity and clustercentric radius cuts of My < -19.8 and 
R < 0.6/?20o) as for the early-type fractions except that the 
early-type selection criteria on bul ge fraction and image s mooth- 
ness were simply replaced by the Poggianti et all d2006l) cut of 
EW([OII]) < — 3A. In order to evaluate the importance of the 
errors on our equivalent widths on our determination of the frac- 
tions of [Oil] emitters, we also computed [Oil] fractions using 
equivalent widths from Brinchmann et al. (2004). The two sets 
of equivalent widths are plotted against one another in Figure [6] 
The agreement between the two sets is excellent, and we con- 
clude that our [Oil] fractions are robust. 

Table|4]gives corrected early-type galaxy fractions and frac- 
tions of [Oil] computed for R < 0.6/? 2 oo for the 158 SDSS clus- 
ters in our local comparison sample. We included only galax- 
ies brighter than Myjim = -19.8 to avoid incompleteness in 
the SDSS spectroscopic sample. This cutoff magnitude corre- 
sponds to the absolute magnitude limits we used for our distant 
EDisCS clusters once passive evolution is taken into account 
(see Section l4~2l i. 



code described in lDe Lucia & Blaizotl d2"007i|pl The MS followed 
2160 3 particles of mass 8.6xlO 8 /i~'M within a comoving box 
of size 500 \C X Mpc on a side with a spatial resolution of 5h~ l 
kpc. Early-type galaxy fractions were computed from these sim- 
ulated galaxy catalogs using the following procedure. Haloes 
were randomly selected at three different redshifts (z = 0, 0.41, 
0.62) so that they were uniformly distributed in log(M2oo)- The 
final halo sample was 100 haloes at z — 0, 94 haloes at z = 0.41 
and 92 haloes at z = 0.62. For each of these haloes, all galax- 
ies in a cubic box 6 Mpc on a side around the central galaxy 
were selected, and a morphological type was assigned to each 
model galaxy by computing the quantity AM = M^uige - M tota i 
(in the rest-frame B-band). Galaxies with AM < 1.0 were con- 
sidered to be "early-type". This is the same criterion as select- 
ing real galaxies with B/Tfousi t 0.40. It is important here 
to note that an early-type galaxy in the simulations was defined 
solely based on this cut in bulge fraction because the simulations 
do not have the resolution required to model internal fine struc- 
tures such as asymmetries. Given that real, early-type galaxies 
were also selected according to image smoothness, one might 
find the early-type fractions of real clusters to be systematically 
lower. For each halo, the fraction of early-type galaxies within 
0.6/?2oo from the BCG was computed using three different pro- 
jections. Furthermore, only galaxies that were within 2 Mpc 
from the BCG along the line of sight were included. The frac- 
tions were computed using only galaxies brighter than -20.5, 
-20.1, and -19.8 in the rest-frame V-band at redshift 0.6, 0.4, 
and 0.0 respectively to match the limits used for the SDSS and 
EDisCS early-type galaxy fractions. A galaxy in the simulation 
was deemed to be star-forming if its star-formation rate in the 
last timestep of its evolution was not equal to zero. 

Figure [7] shows the resulting model early-type fractions as 
a function of cluster velocity dispersion, redshift and fraction 
of star-forming galaxies for the MS haloes. At a given red- 
shift, there is no dependence of the early-type fraction on cluster 
velocity dispersion, but the scatter symmetrically increases to- 
wards both lower and higher fractions leading to a "wedge-like" 
distribution towards lower cluster cr's. The early-type fractions 
of both low and high-mass clusters increase with decreasing red- 
shift from -0.70 at z = 0.65 to -0.85 at z = 0. The early-type 
fractions of massive clusters are anticorrelated with the fractions 
of star-forming galaxies: clusters at z = have higher early- 
type fractions but lower fractions of star-forming galaxies. Note 
that the trends in F igure [7] do not agree with those shown in 
Diaferio et al. (2001) although the assumptions made about mor- 
phological transformations are very similar in the two models. In 
particular, the MS shows little trend of early-type fraction with 
cluster velocity dispersion but a substantial trend with redshift, 
while Diaferio et al. found the opposite. This is likely a result of 
the poorer mass resolution, poorer statistics and cruder dynami- 
cal modelling of the earlier paper. 



5. Results 



4.5. Theoretical Models 

Numerical simulations of dark matter haloes populated with 
galaxies using semi-analytical models greatly help in the inter- 
pretation of obse rvational results. We use here the Millennium 
Simulation (MS; Springelet al., 2005), and the semi-analytical 



We use here our VLT/FORS2 early-type fractions for all EDisCS 
clusters for the sake of uniformity. 



2 Simulated galaxy catalogs used here are publicly available at 
http://www.mpa-garching.mpg.de/millennium/ 
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Fig. 7. Early-type galaxies in Millennium Simulation dark matter haloes Top, left-hand panel: Early -type galaxy fraction within 
0.6^200 versus cluster velocity dispersion at three different redshifts. Top, right-hand panel: Early-type galaxy fraction within 
0.6/? 200 versus age of the universe. Blue and red points are clusters with velocity dispersions below and above 600 km/s respectively. 
Lower, left-hand panel: Early-type galaxy fraction within 0.6/?200 versus fraction of star-forming galaxies in clusters with cr < 600 
km/s. Blue points show haloes selected at redshift zero, and all the other haloes are in red. Lower, right-hand panel: Early-type 
galaxy fraction within 0.6/?20o versus fraction of star-forming galaxies in clusters with cr > 600 km/s. 



5.1. Early-Type Galaxy Fractions versus Cluster Velocity 
Dispersion and Redshift 

Figure [8] shows early-type galaxy fractions versus velocity dis- 
persion for the SDSS and EDisCS clusters. The early-type 
galaxy fractions of both cluster samples exhibit no clear trend 
as a function of cr. Table [5] gives Spearman rank test results for 
the SDSS sample and different EDisCS subsamples. The only 
significant correlation between early-type fraction and velocity 
dispersion is found in the high-z EDisCS clusters. It only has a 
2.5% chance of being due to r andon sampl i ng. S uch a positive 
correlation was also reported in lDesai et alj d2007l) for the same 
cluster subsample, but it disappears when the full EDisCS sam- 
ple is considered. The lack of a significant correlation agrees 
well with the results for the Millennium Simulation in the top 
left-hand panel of Figure Q but disagrees with the earlier the- 



oretical results of Diaferio et al. (2001) which showed a trend 
between f et and cr. A visual inspection of Figure [8] confirms 
the statistical test results. The mid-z EDisCS clusters do not 
show any correlation with cr in contrast to the high-z clusters. 
In particular, two mid-z EDisCS clusters (CL1 1 19.3-1 130 and 
CL1420. 3-1236) with cr ~ 200 km/s have early-type fractions 
similar or higher (f e , ~ 0.5-0.8) than the most massive clusters 
in our sample. Interestin gly, the same two clusters were found 
bv lPoggianti et al.l d2006) to be the most outstanding outliers in 
the [Oil] fraction - cr relation in the sense that they have a low 
fraction of [Oil] emitters for their mass. This is consistent with 
what we observe here given that early-type galaxies typically 
have lower [Oil] emission fluxes. 

Figure [8] does show that there is a marked difference in the 
morphological content of the EDisCS and SDSS clusters. All 
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Table 5. Spearman Rank Tests Results for Early-Type Fraction 
versus Cluster Velocity Dispersion 



Cluster Sample 


Na 


R s 


p- value 


(1) 


(2) 


(3) 


(4) 


SDSS 


158 


-0.05 


0.51 


EDisCS all 


18 


0.18 


0.47 


EDisCS mid-z 


9 


-0.11 


0.78 


EDisCS high-z 


9 


0.73 


0.025 



Table 6. Two-sample Kolmogorov-Smirnov Test Probabilities 
for Early-Type Fraction versus Cluster Velocity Dispersion 





SDSS 


Cluster Sample 2 
SDSS EDisCS 


EDisCS 


Cluster Sample 1 


(All) 


{a > 600) (a < 600) 


(o- > 600) 


SDSS (a < 600) 




0.628 0.190 




SDSS (cr > 600) 






0.472 


EDisCS (a < 600) 






0.250 


EDisCS (All) 


0.506 







EDisCS f et values (with the exception of one cluster) are be- 
low 0.6, but half of the SDSS clusters are above this value. The 
population of early-type galaxies has thus increased significantly 
in half of the clusters of all velocity dispersions. An increase 
in early-type fraction with decreasing redshift may already be 
visible when one compares mid-z and high-z EDisCS clusters. 
Mid-z clusters around cr = 500 km/s have f et 0.5 whereas the 
high-z clusters have f et 0.4. This would represent a ~25% in- 
crease over a time interval of 2 Gyrs. As shown in Figure |7l the 
early-type fractions of clusters in the Millenium Simulation also 
increase with decreasing redshift in clusters of all velocity dis- 
persions, but there is a lack of simulated clusters with f e , < 0.5 
compared with the SDSS-C4 clusters. The scatter in the f et val- 
ues of simulated clusters is also smaller than in those of real 
clusters. For simulated clusters at z = with cr > 600 km/s, 
o~(fet) = 0.06 compared to o-(f et ) = 0.21 for SDSS clusters over 
the same range of velocity dispersions. Given that the mean er- 
ror on the SDSS f et values is 0.12, the intrinsic scatter would be 
0.17. This intrinsic scatter is still almost three times the scatter 
in the simulated clusters. 

Figure [9] shows SDSS and EDisCS early-type fractions as a 
function of the age of the universe (i.e., redshift). The clusters 
have been divided into two subgroups based on their velocity 
dispersions. The early-type fractions of massive {cr > 600 km/s) 
EDisCS clusters (right panel) are in very good ag reement with 
the ones in the compilation of Ivan Dokkum et al] d200ll) which 
also have velocity disper sions greater than 600 km /s. The clus- 
ters at low redshift in the van Dokkum et al] d2001l) compilation 
suggest that there are no local clusters with low early-type frac- 
tions and hence that all clusters have uniformly increased their 
early-type fraction from z ~ 0. to the present day. However, our 
SDSS cluster sample shows that this simple picture is not en- 
tirely true. While half of the SDSS clusters have higher early- 
type fractions than clusters at high redshift, the other half have 
early-type fractions equal or even lower than the EDisCS clus- 
ters. The same holds true for the low mass clusters (left-hand 
panel). The scatter in f et (< 0.1) in high-mass EDisCS clusters 
does appear to be considerably less that the scatter seen in low- 
mass clusters. 

The lack of a clear trend in early-type fraction with redshift 
in the right-hand panel of Figure [9] is in disagreement with the 
Millennium Simulation prediction in the top right-hand panel of 



Table 7. Spearman Rank Tests Results for Early-Type Fraction 
versus Fraction of [Oil] Emitters 



Cluster Sample 


N d 


R s 


p- value 


(1) 


(2) 


(3) 


(4) 


SDSS (All) 


158 


-0.63 


4.33x10-'" 


SDSS (cr < 600) 


108 


-0.57 


1.50x10-'° 


SDSS (cr > 600) 


50 


-0.77 


6.04x10-" 


EDisCS (All) 


18 


-0.74 


0.00043 


EDisCS (cr < 600) 


11 


-0.78 


0.0043 


EDisCS (cr > 600) 


7 


-0.77 


0.0438 



Table 8. Two-sample Kolmogorov-Smirnov Test Probabilities 
for [Oil] Emitter Fraction versus Cluster Velocity Dispersion 





SDSS 


Cluster Sample 2 
SDSS EDisCS 


EDisCS 


Cluster Sample 1 


(All) 


(tr > 600) (o- < 600) 


(cr > 600) 


SDSS (cr < 600) 




0.090 0.005 




SDSS (cr > 600) 






0.046 


EDisCS (cr < 600) 






0.761 


EDisCS (All) 


0.026 







Figure [7] There is a clear deficit of clusters with low early-type 
fraction at low redshift in the Millenium Simulation compared 
with our SDSS sample. 

5.2. Early-Type Galaxy Fractions versus Fractions of [Oil] 
Emitters 

The link between star formation and morphological transforma- 
tion and its evolution as a function of redshift provides more 
clues on the processes driving galaxy morphology in local and 
distant clusters. The fractions of galaxies with [OH ] emission 
in the EDisCS clusters were computed as in Pog gianti et al] 
(2006) using the same absolute magnitude limits and the same 
prescriptions for correcting magnitude and geometric incomplet- 
ness, but the clustercentric radius cut was changed to match 
the one used for the early-type fractions in this paper (R et < 
0.6i?20o)' The two datasets are therefore directly comparable. 
Figure[10]shows f et versus f[oii] with our local and distant sam- 
ples again divided according to velocity dispersion. Table|7]gives 
Spearman test results between f et and f[pin- There is a strong 
correlation between f et versus f[oii] in both SDSS and EDisCS 
cluster samples irrespective of cluster velocity dispersion. The 
EDisCS clusters lie within the envelopes defined by the SDSS 
clusters. There is no offset between the zeropoints of the corre- 
lations at low and hig h redshift. However, as demonstrated by 
Poggianti et alj d2006l) . the star formation activity (parametrized 
by f[oii]) has decreased in all environments from z ~ 0.75 to 
z ~ 0.08. This is confirmed by the K-S test results in Table [8] 
The probabilities that the EDisCS and SDSS clusters are drawn 
from the same parent f[oir\ distribution are only 0.026, 0.005 
and 0.046 for the whole samples, low cr and high cr subsamples 
respectively. 

The f et versus foil] values for clusters from the Millenium 
Simulation (Figure 01 are quite different from the observations. 
Low cr MS clusters at low and high redshifts are confined to high 
f et and fon values with no apparent correlation. There is only a 
handful of clusters with low values for both f et and fon- The 
high cr MS clusters are found in a very limited range of f e , and 
fon values (0.35 < f n < 0.75, 0.6 < f et < 0.85). 



16 



Simard et al.: Evolution of the Early-Type Galaxy Fraction in Clusters 



- 1 1 1 r- 



~1 ' ' ' ' 1 1- 



• . • • •••• • 



••• • • 



_i i i i 



_i i i i 



N 
VI 
OS 

o 



500 
a (km/s) 



1000 



—\ 1 1 1 - 



I 1 1 1— 



_i I I 1_ 



_i I I I 



500 
a (km/s) 



_i_ 



1000 



Fig. 8. Early-type galaxy fraction within O.6R200 versus velocity dispersion for SDSS and EDisCS clusters. Both samples have been 
matched in velocity dispersion. Left panel: SDSS clusters. Only typical error bars are shown in the lower right-hand corner for 
clarity. Right panel: Filled and open circles are mid-z and high-z EDisCS clusters respectively. Errors bars shown in both panels are 
lcr errors. Our VLT/FORS2 early-type fractions are used here for all EDisCS clusters for the sake of uniformity. 
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Fig. 9. Early-type galaxy fraction versus age of the universe (i.e., redshift) for clusters with cr < 600 km/s (left panel) and clusters 
with cr > 600 km/s (right panels). SDSS and EDisCS clusters are blue an d red respectively, a nd both samples have been matched in 
velocity dispersion. Clusters shown in black are from the compilation of Ivan Dokkum et aD (1200 ll) in which open and solid points 
have X-ray luminosities below and over 10 44 5 ergs s _1 respectively. Our VLT/FORS2 early-type fractions are used here for all 
EDisCS clusters for the sake of uniformity. 



6. Discussion 

In order to fully understand possible evolutionary trends ob- 
served here, it is important to determine how cluster velocity 
dispersion changes with redshift as a result of the hierarchical 
growth of structures. Are we looking at similar clusters when 
we focus on the same range of velocity disp ersions in the SDSS 
and EDisCS clusters? Poggi anti et alj ([2006) looked at the mean 
change in cr between z — and z = 0.76 using a sample of 90 
haloes from the Millennium Simulation uniformly distributed in 
log(mass) between 5 x 10 12 and 5 x 10 15 M Q . Their Figure 8 
shows how cr evolves over that redshift interval. For example, 
a z = cluster with cr = 900 km/s would typically have cr ~ 
750 km/s at z - 0.76. This evolution is not sufficient to introduce 
biases in our analysis here. Indeed, selecting clusters with cr > 
600 km/s, say, at either z — or z — 0.76 would keep nearly all 
the same clusters. Measured velocity dispersions may exhibit a 



large scatter with respect to the true halo mass particularly for 
low-mass clusters. The velocity dispersions for the SDSS and 
EDisCS clusters were calculated in a very similar way in order 
to minimize any biases. Velocity dispersions calculated from a 
small number of cluster members may be overestimates of the 
true cluster mass. TableQ]lists 1 103.7-1245b as the cluster with 
the lowest number of members (N — 1 1). In order to check the 
robustness of our results, we re-ran our analyses by excluding 
SDSS clusters in Table|4]with N < 10 for which velocity disper- 
sions may be less reliable and found that our results remained 
unchanged. 

Pogg ianti et alj d2006l) proposed a scenario in which two 
channels are responsibl e for the production of passive galaxies 
in clusters, and others (Fab er et all 120071: iBrown et al.L [2007) 
have proposed a similar scenario for the migration of galaxies 
from the "blue cloud" to the red sequence. "Primordial passive 
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Fig. 10. Early-type galaxy fraction versus [Oil] emitter fraction for clusters with cr < 600 km/s (left panel) and clusters with 
cr > 600 km/s (right panel). SDSS and EDisCS clusters are shown in blue and red respectively, and both samples have been 
matched in velocity dispersion. Only typical error bars are shown for the SDSS clusters in the lower right-hand corner for clarity. 
Our VLT/FORS2 early-type fractions are used here for all EDisCS clusters for the sake of uniformity. 



galaxies" are composed of galaxies whose stars all formed at 
very high redshift (z > 2) over a short timescale. These galax- 
ies have been observed in clusters up and beyond z = 1, and 
they largely comprise luminous ellipticals. "Quenched passive 
galaxies" have had a more extended period of star formation 
activity, and their star formation has been quenched after their 
infall into dense cluster environments. These quenched passive 
galaxies would then suffer the effects of cluster processes such 
as ram pressure stripping, harassment, strangulation and merg- 
ers to become SO and earlier type galaxies. A key point of this 
scenario is that processes affecting morphology and star forma- 
tion activity operate on d ifferent timescales as shown rece ntly for 
the EDisCS sample by ISanchez-Blazquez et al.l (120091) . There 
is good evidence that star formation is quenched in galaxies 
over timescal es of 1-3 Gyr after they hav e entered the cluster 
environment (PoggiantreTafl, [1999[ [2006) whereas morpholog- 
ical transformati on through mergers and harassment can take 
longer (~ 5 Gyr, iMoore et al.L 119981) . The best example of this 
is the fact that the vast majority of post-starburst galaxies in 
distant clusters, those that have had their star formation ac- 
tivity ter minated during t he las t Gyr, still retain a spiral mor- 
phology (Pogg ianti et all 1 19991) . Such a two-channel scenario 
would naturally explain observations indicating that the ellipti- 
cal galaxy fraction actually remains constant wit h redshift while 
the SO fraction rises with decreasing reds hift dDressler et all 
ll997HFasano et alll2000tlDesai et aUl2007l) . Unfortunately, the 
VLT/FORS2 images do not have sufficient spatial resolution to 
disentangle E and SO galaxies as mentioned in Section PTTl to de- 
termine the exact contribution from each channel. We can there- 
fore only study the overall production of early-type galaxies, but 
it should exhibit different behaviors with cluster global proper- 
ties depending on the process(es) dominating it. Given our quan- 
titative definition of an early-type galaxy based on bulge fraction 



and image smoothness, there are essentially two ways to trans- 
form late-type galaxies into early-type ones: 1) processes such as 
collisions and harassment that can fundamentally alter the struc- 
ture of a galaxy by forming bulges and/or destroying disks and 
2) quenching processes that can extinguish star forming regions 
responsible for some of the galaxy image asymmetries and also 
cause a fading of t he disks. 

Applying the iPoggianti et all (120061) scenario to our re- 
sults, the "threshold" in f et values in our high redshift clusters 
(Figures [8] and |9]l could be explained by a population of pri- 
mordial passive galaxies that formed at even higher redshifts. 
Most of our high redshift clusters have early-type fractions in 
the range 0.3-0.6 with no correlation with cluster velocity dis- 
persion. Are these early-type fractions indeed consistent with a 
po pulations of primordia l passive galaxies? Calculations done 
in IPoggianti et al. (2006) show that the fraction of galaxies at 
Z = 0.6 that were present in haloes with masses greater than 
3xl0 12 M at z = 2.5 is 0.4+0.2. These primordial passive galax- 
ies can therefore account for at least 2/3 (if not all) of the early- 
type populations in high redshift clusters, and their high forma- 
tion redshift would explain the lack of dependence of f et on clus- 
ter velocity dispersion. 

One of our main results is that the early-type fractions of 
galaxy clusters increase from z = 0.6 - 0.8 to z ~ 0.08 in 
clusters of all velocity dispersions. What kind of morpholog- 
ical transformation process(es) can lead to such an evolution? 
Collisions and harassment both depend on galaxy-galaxy inter- 
actions and the time a galaxy has spent within the cluster en- 
vironment. Cluster velocity dispersion influences the number of 
interactions and their duration. Higher velocity dispersions in 
more massive clusters yield more interactions per unit time N 
but with shorter durations Af in a given time interval. One might 
therefore expect to see a peak in early-type type fraction at the 
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cluster velocity dispersion where the product NAf is maximized. 
No such peak is seen in our clusters. Ram-pressure stri pping 
is expected to go as (nrcMV 2 '^)/Mrep dGaetz et all Il987l) with 

kicm, Vg a i and M rep being the density of the ICM, the veloc- 
ity of the galaxies within the ICM and the rate at which galax- 
ies can replenish their gas respectively. The fraction of passive 
galaxies should therefore be a relatively strong function of clus- 
ter velocity dispersion if quenching by ram pressure stripping is 
the dominant process. The number of post-starburst galaxies in 
EDisCS clusters does c orrelate with cluster velocity dispersion 
(Poggianti et all l2009al) . but the uniform increase in early-type 
fractions at all cluster velocity dispersions observed going from 
EDisCS to SDSS clusters is not consistent with the intracluster 
medium being the main cause of the changes in cluster morpho- 
logical content. 

Even though the early-type and [Oil] emitter fractions in 
EDisCS and SDSS clusters show no corr elation with cluster ve- 
locity dispersion (Pogg ianti et all [20061 and this work), there 
is a very strong correlation between f et and fan. This corre- 
lation is seen at both low and high cluster masses as well as 
at both low and high redshifts. Morphology and star formation 
therefore appear to be closely linked with one another over a 
wide range of environments and times. However, different struc- 
tural transformation and quenching processes are thought to 
opera te over different timescales (e.g., San chez-Blazquez et all 
2009). Timescales range from 1-2 Gyr (based on typical clus- 
ter crossing times) for truncating star formation to 3-5 Gyr for 
totally extinguishing star formation in new l y acc reted galaxies 
dPoggianti et all l2006t iTonnesen & Brvanl l2009l) . Looking at 
the evolution of EDisCS cluste r red-sequence galaxies over 2 
Gyr (from z = 0.75 to z = 0.451. ISanchez-Blazquez eTail (2009) 
found that morphological transformation and quenching of star 
formation indeed appeared to not be simultaneous. As noted 
in Section 15.11 the early-type fractions of mid-z EDisCS clus- 
ters may be ~25% higher than the ones of high-z clusters. This 
change would therefore have taken place over a 2 Gyr interval 
in our adopted cosmology. However, the time baseline here be- 
tween SDSS and EDisCS clusters is almost 6 Gyr, and, unfor- 
tunately, this is ample time to erase any difference arising from 
different timescales in the link between morphology and star for- 
mation. 

The lack of dependence of morphology and star forma- 
tion on global cluster properties such as velocity dispersion 
raises the question of whether changes in galaxy properties are 
driven by more local effects or whether they occur out s ide of 



the cluster environment. Recent w ork (Poggianti et al 
iPark & ChoiL l2009t iBamford et all l2009t lEllison et al" 



have re-emphasized the strong link between galaxy properties 
and local galaxy density rather than cluster membership. Galaxy 
properties are seen to change at densities around 15-40 galax- 
ies Mpc~ 2 or projected separations of 20-30 /r 1 kpc. Others 
(e.g., iKautsch et all 120081: IWilman et all 120091) have suggested 
that the galaxy group environment might be more conducive 
to galaxy transformation. Our observed evolution in early-type 
fraction as a function of redshift and the strong correlation be- 
tween morphology and star formation at all cluster masses would 
support the idea that cluster membership is of lesser importance 
than other variables such as local density in determining galaxy 
properties. 

The properties of simulated clusters from the Millenium 
Simulation compare well with those of EDisCS and SDSS clus- 
ters. Their early-type fractions also show no dependence with 
cluster velocity dispersion in contrast to previous theoretical 



work (e.g. iDiaferio et all 1200 lb but in agreement with observa- 
tions. However, there is a definite lack of MS clusters with low 
early-type fractions at z — compared to the SDSS sample. It 
is important here to note that an early-type galaxy in the sim- 
ulations was defined solely based on its bulge fraction because 
the simulations do not have the resolution required to model in- 
ternal fine structures such as asymmetries. Given that real, early- 
type galaxies were also selected according to image smoothness, 
one would expect the early-type fractions of real clusters to be 
systematically lower. However, half of the SDSS clusters have 
low early-type fractions not seen in the simulations at z — 0, and 
such a large discrepancy could only be explained by a signifi- 
cant population of real bulge-dominated galaxies with relatively 
large asymmetries. It is more likely that bulge formation in the 
simulations may be too efficient. The scatter in f et values for the 
simulated clusters with cr > 600 km/s is also nearly three times 
smaller than observed in the real clusters (Section 15. U which 
may indicate that the models may not include the right mixture 
of evolutionary processes at work on real galaxies. High-mass 
simulated clusters show a correlation between early-type frac- 
tion and star-forming fraction (albeit over narrower ranges than 
observed), but the correlation is not seen in the low-mass sim- 
ulated clusters. This may be understood by high mass clusters 
having been formed long enough for evolutionary processes to 
have had enough time to act on galaxies to modify their proper- 
ties whereas this is not necessarily the case for low-mass clus- 
ters. The fact that the correlation is observed in both low- and 
high-mass real clusters may be an indication that processes giv- 
ing rise to the correlation may be more efficient (or altogether 
different) than modelled. It is also important to keep in mind 
here that the properties of a galaxy in these models are essen- 
tially driven by the mass of its parent halo. 



7. Summary 

We have presented quantitative morphologies measured from 
PSF-convolved, 2D bulge+disk decompositions of cluster and 
field galaxies on deep VLT/FORS2 images of eighteen, 
optically-selected galaxy clusters at 0.45 < z < 0.80 observed as 
part of the ESO Distant Cluster Survey. The morphological con- 
tent of these clusters was characterized by the early-type fraction 
within a clustercentric radius of 0.67? 200, and early-type galaxies 
were selected based on bulge fraction and image smoothness. We 
showed a very good agreement between quantitative and visual 
galaxy classifications. We used a set of 158 clusters extracted 
from the Sloan Digital Sky Survey matched in velocity disper- 
sion to our EDisCS sample and analyzed exactly in the same 
way to provide a robust comparison baseline and to control sys- 
tematics. We studied trends in early-type fraction as a function 
of cluster mass and redshift. We also explored the link between 
morphology and star formation by comparing early-type frac- 
tions to the fractions of [Oil] emitters in our clusters. Our main 
results are: 

1. The early-type fractions of the SDSS and EDisCS clusters 
exhibit no clear trend as a function of cluster velocity dispersion. 

2. Mid-z EDisCS clusters around cr = 500 km/s have f e , 0.5 
whereas high-z EDisCS clusters have f et 0.4. This represents 
a ~25% increase over a time interval of 2 Gyrs. 

3. There is a marked difference in the morphological content 
of the EDisCS and SDSS samples. None of the EDisCS clus- 
ters have an early-type fraction greater than 0.6 whereas half of 
the SDSS clusters lie above this value. This difference is seen in 
clusters of all velocity dispersions (i.e., masses). 
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4. There is a strong and clear correlation between morphol- 
ogy and star formation activity in the sense that decreasing frac- 
tions of [Oil] emitters are tracked by increasing early-type frac- 
tions. This correlation holds in both low and high cluster masses 
as well as at both low and high redshift. 

5. The early-type fractions of clusters drawn from the 
Millennium Simulation (Springel et al., 2005) using the galaxy 
formation model of |Pe Lucia & Blaizotl (|2007) also show no 
clear dependence on cluster velocity dispersion. However, at z 
= 0, they are not enough simulated clusters with low early-type 
fractions compared to the SDSS cluster sample. While high- 
mass simulated clusters show a correlation between early-type 
fraction and star-forming fraction (albeit over narrower ranges 
than observed), this correlation is not seen in the low-mass sim- 
ulated clusters in contrast to the real ones. 

Our results pose an interesting challenge to structural trans- 
formation and star formation quenching processes that strongly 
depend on the global cluster environment (e.g., a dense ICM) 
and suggest that cluster membership may be of lesser impor- 
tance than other variables in determining galaxy properties. 
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3041 


0.0776 


15 


555+72 


1.32 


0.283 


0.45+0.15 


0.46+0.15 


0.45+0.15 


0.45+0.15 


2016 


0.0449 


135 


552+34 


1.34 


0.283 


0.42+0.09 


0.41+0.09 


0.33+0.08 


0.31+0.08 


3474 


0.0519 


19 


547+93 


1.32 


0.275 


0.32+0.15 


0.27+0.16 


0.44+0.16 


0.48+0.17 


3505 


0.0711 


18 


543+96 


1.30 


0.266 


0.55+0.15 


0.55+0.15 


0.45+0.15 


0.50+0.15 


3237 


0.0456 


30 


542+66 


1.31 


0.267 


0.85+0.10 


0.90+0.10 


0.32+0.13 


0.30+0.14 


3249 


0.0740 


19 


541+98 


1.29 


0.263 


0.60+0.13 


0.63+0.13 


0.47+0.13 


0.45+0.13 


3117 


0.0623 


49 


536+52 


1.29 


0.257 


0.72+0.10 


0.75+0.10 


0.13+0.07 


0.08+0.07 


3003 


0.0597 


59 


535+56 


1.29 


0.255 


0.61+0.09 


0.62+0.09 


0.35+0.09 


0.35+0.09 


3121 


0.0627 


29 


529+71 


1.27 


0.247 


0.63+0.12 


0.62+0.12 


0.30+0.11 


0.31+0.12 


3167 


0.0732 


19 


526+82 


1.26 


0.240 


0.70+0.13 


0.71+0.13 


0.38+0.13 


0.37+0.13 


2058 


0.0734 


21 


526+99 


1.26 


0.241 


0.54+0.14 


0.52+0.14 


0.38+0.13 


0.39+0.14 


3011 


0.0820 


26 


524+81 


1.25 


0.237 


0.66+0.11 


0.66+0.11 


0.24+0.10 


0.23+0.10 


3015 


0.0689 


34 


523+67 


1.25 


0.238 


0.59+0.10 


0.59+0.10 


0.36+0.10 


0.35+0.10 


2178 


0.0526 


14 


522+84 


1.26 


0.238 


0.26+0.17 


0.21+0.17 


0.58+0.19 


0.60+0.19 


1047 


0.0829 


22 


521+87 


1.24 


0.233 


0.92+0.09 


1.00+0.09 


0.08+0.09 


0.00+0.09 


1120 


0.0756 


13 


516+109 


1.23 


0.227 


0.44+0.16 


0.41±0.17 


0.44+0.16 


0.44+0.17 


1050 


0.0717 


14 


514+88 


1.23 


0.225 


0.68+0.13 


0.65+0.14 


0.41+0.14 


0.44+0.14 


3583 


0.0673 


29 


512+80 


1.23 


0.223 


0.64+0.10 


0.63+0.10 


0.31+0.10 


0.31+0.10 


3196 


0.0603 


17 


512+92 


1.23 


0.224 


0.29+0.14 


0.26+0.14 


0.50+0.16 


0.51+0.16 


2005 


0.0433 


48 


506+47 


1.23 


0.218 


0.84+0.11 


0.91+0.11 


0.16+0.11 


0.11+0.11 


1223 


0.0600 


18 


506+50 


1.22 


0.216 


0.50+0.34 


0.74+0.34 


0.50+0.34 


0.41+0.34 


2020 


0.0740 


32 


505+83 


1.21 


0.213 


0.68+0.15 


0.71+0.15 


0.08+0.09 


0.00+0.09 


1179 


0.0775 


16 


505+94 


1.20 


0.212 


0.50+0.14 


0.45±0.14 


0.50+0.14 


0.56+0.14 


3539 


0.0542 


12 


502+122 


1.21 


0.212 


0.21+0.20 


0.00+0.20 


0.50+0.25 


0.56+0.25 


3080 


0.0616 


17 


501+110 


1.20 


0.209 


0.68+0.15 


0.74+0.16 


0.32+0.15 


0.25+0.16 


2111 


0.0778 


8 


500+172 


1.19 


0.207 


0.50+0.34 


0.77+0.34 


0.50+0.34 


0.00+0.34 


1149 


0.0523 


33 


500+62 


1.21 


0.209 


0.45+0.15 


0.46+0.15 


0.55+0.15 


0.53+0.15 


1404 


0.0789 


14 


498+155 


1.19 


0.204 


0.44+0.16 


0.35+0.16 


0.56+0.16 


0.66+0.16 


1144 


0.0751 


20 


497+71 


1.19 


0.203 


0.71+0.14 


0.71+0.14 


0.29+0.14 


0.29+0.14 


1147 


0.0809 


17 


493+95 


1.17 


0.197 


0.68+0.13 


0.69+0.13 


0.15+0.10 


0.14+0.10 


3019 


0.0682 


32 


492+67 


1.18 


0.198 


0.74+0.11 


0.78+0.11 


0.53+0.13 


0.54+0.13 


1178 


0.0543 


32 


490+41 


1.18 


0.197 


0.54+0.14 


0.51+0.14 


0.54+0.14 


0.60+0.13 


3066 


0.0767 


30 


482+49 


1.15 


0.185 


0.61+0.10 


0.62±0.10 


0.30+0.09 


0.24+0.09 


2027 


0.0790 


34 


481+62 


1.14 


0.184 


0.50+0.34 


0.54+0.34 


0.50+0.34 


0.18+0.34 


3344 


0.0738 


8 


479+63 


1.14 


0.182 


0.80+0.13 


0.78+0.16 


0.32+0.15 


0.35+0.16 


3083 


0.0760 


19 


470+64 


1.12 


0.172 


0.50+0.13 


0.53±0.13 


0.50+0.13 


0.48+0.13 


1087 


0.0782 


15 


465+180 


1.11 


0.166 


0.37+0.12 


0.37+0.12 


0.63+0.12 


0.60+0.12 


1005 


0.0809 


21 


458+50 


1.09 


0.159 


0.50+0.20 


0.50+0.20 


0.31+0.19 


0.25+0.19 


1360 


0.0746 


8 


448+103 


1.07 


0.149 


0.23+0.15 


0.18+0.15 


0.77+0.15 


0.84+0.15 


3120 


0.0731 


19 


445+80 


1.06 


0.146 


0.70+0.13 


0.70±0.13 


0.38+0.13 


0.38+0.13 


3248 


0.0446 


20 


442+90 


1.07 


0.145 


0.26+0.17 


0.27±0.17 


0.58+0.19 


0.56+0.19 


3069 


0.0555 


20 


441+47 


1.06 


0.144 


0.76+0.12 


0.82+0.12 


0.41+0.14 


0.39+0.14 


2121 


0.0817 


13 


441+86 


1.05 


0.141 


0.42+0.19 


0.37+0.19 


0.58+0.19 


0.63+0.19 


2141 


0.0779 


16 


422+89 


1.00 


0.124 


0.80+0.13 


0.87+0.13 


0.20+0.13 


0.11+0.13 


3506 


0.0635 


11 


421+87 


1.01 


0.124 


0.26+0.17 


0.22±0.17 


0.74+0.17 


0.79+0.17 


1368 


0.0520 


10 


417+91 


1.01 


0.121 


0.61+0.22 


0.72+0.22 


0.39+0.22 


0.28+0.22 


2052 


0.0772 


12 


415+72 


0.99 


0.118 


0.50+0.14 


0.49±0.14 


0.15+0.10 


0.10+0.10 


2000 


0.0699 


29 


411+53 


0.98 


0.115 


0.48+0.10 


0.48+0.10 


0.43+0.10 


0.43+0.10 


1009 


0.0746 


27 


404+41 


0.97 


0.110 


0.53+0.11 


0.52+0.11 


0.31+0.11 


0.33+0.11 


3140 




Z4 


4Uj±4j 




A 1 1 A 

U.l 1U 


U.jj+U.Ij 


0.53+0.15 


0.07+0.08 


<J.<J<J+<J.<JO 


1043 


0.0743 


30 


403+60 


0.96 


0.109 


0.74+0.11 


0.79+0.11 


0.26+0.11 


0.22+0.11 


1106 


0.0402 


18 


402+112 


0.98 


0.109 


0.39+0.15 


0.34+0.15 


0.50+0.16 


0.44+0.16 


3572 


0.0587 


11 


393+183 


0.95 


0.102 


0.79+0.20 


1.00+0.20 


0.21+0.20 


0.00+0.20 


3144 


0.0805 


16 


386+57 


0.92 


0.095 


0.94+0.06 


1.00+0.06 


0.14+0.09 


0.08+0.10 


2010 


0.0645 


18 


384+55 


0.92 


0.095 


0.50+0.18 


0.46+0.18 


0.36+0.17 


0.36+0.17 


2211 


0.0422 


22 


372+38 


0.90 


0.087 


0.68+0.15 


0.74+0.16 


0.32+0.15 


0.31+0.16 
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SDSS C4 


z 


N 


cr v 


^200 


M d 




R < 0.6X 2 oo 




ID 






(km/s) 


(Mpc) 


(10 15 M o ) 


fet,raw 


fet,corr 


flOm.raw 


f[OII],corr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1341 


0.0700 


10 


370+59 


0.88 


0.084 


0.58±0.19 


0.56±0.19 


0.26±0.17 


0.21±0.17 


3388 


0.0722 


12 


357±109 


0.85 


0.075 


0.61 ±0.22 


0.66±0.22 


0.61 ±0.22 


0.63±0.22 


3007 


0.0685 


32 


356+50 


0.85 


0.075 


0.73±0.09 


0.74±0.09 


0.17±0.08 


0.13±0.08 


2140 


0.0810 


5 


355+216 


0.84 


0.074 


0.13±0.14 


0.00±0.14 


0.87±0.14 


1.00±0.14 


3182 


0.0614 


18 


352+84 


0.85 


0.073 


0.44±0.16 


0.42±0.17 


0.20±0.13 


0.16±0.13 


1006 


0.0477 


24 


340+54 


0.82 


0.066 


0.77±0.15 


0.85±0.15 


0.23±0.15 


0.15±0.15 


2120 


0.0518 


12 


338+101 


0.82 


0.065 


0.58±0.19 


0.58±0.19 


0.42±0.19 


0.41±0.19 


3341 


0.0723 


14 


331+97 


0.79 


0.060 


0.56±0.16 


0.55±0.17 


0.32±0.15 


0.30±0.16 


3272 


0.0692 


9 


329+54 


0.79 


0.059 


0.69±0.19 


0.72±0.19 


0.50±0.20 


0.55±0.20 


2156 


0.0668 


17 


328+50 


0.79 


0.059 


0.82±0.12 


0.89±0.12 


0.18±0.12 


0.11±0.12 


1015 


0.0792 


10 


327+66 


0.78 


0.058 


0.79±0.20 


1.00±0.20 


0.50±0.25 


0.57±0.25 


2122 


0.0826 


7 


315+61 


0.75 


0.052 


0.42±0.19 


0.31±0.19 


0.42±0.19 


0.31±0.19 


3340 


0.0611 


8 


314+60 


0.75 


0.051 


0.79±0.20 


1.00±0.20 


0.21 ±0.20 


0.00±0.20 


3247 


0.0571 


14 


307+97 


0.74 


0.048 


0.29±0.26 


0.00±0.26 


0.71 ±0.26 


1.00±0.26 


3061 


0.0480 


21 


300+37 


0.72 


0.045 


0.58±0.19 


0.58±0.19 


0.42±0.19 


0.37±0.19 


1040 


0.0539 


16 


299±38 


0.72 


0.045 


0.39±0.22 


0.40±0.22 


0.39±0.22 


0.36±0.22 


3267 


0.0710 


5 


298+147 


0.71 


0.044 


0.31±0.19 


0.20±0.19 


0.50±0.20 


0.50±0.20 


3328 


0.0412 


12 


297 ±40 


0.72 


0.044 


0.50±0.25 


0.54±0.25 


0.50±0.25 


0.43±0.25 


2070 


0.0775 


12 


295+52 


0.70 


0.043 


0.56±0.16 


0.57±0.17 


0.44±0.16 


0.42±0.17 


3169 


0.0756 


7 


291 ±146 


0.70 


0.041 


0.32±0.15 


0.30±0.16 


0.80±0.13 


0.91 ±0.1 3 


3262 


0.0452 


32 


282±48 


0.68 


0.038 


0.61 ±0.22 


0.67±0.22 


0.61 ±0.22 


0.74±0.22 


1297 


0.0451 


10 


271 ±85 


0.66 


0.033 


0.39±0.22 


0.27±0.22 


0.39±0.22 


0.33±0.22 


3147 


0.0534 


1 1 


265±80 


0.64 


0.031 


0.69±0.19 


0.82±0.19 


0.13±0.14 


0.00±0.14 


1128 


0.0619 


12 


257±64 


0.62 


0.028 


0.64±0.17 


0.66±0.17 


0.36±0.17 


0.35±0.17 


3062 


0.0615 


12 


256±44 


0.62 


0.028 


0.56±0.16 


0.55±0.17 


0.44±0.16 


0.45±0.17 


3318 


0.0419 


7 


252±56 


0.61 


0.027 


0.50±0.20 


0.48±0.21 


0.50±0.20 


0.51 ±0.20 


1171 


0.0556 


15 


242 ±48 


0.58 


0.024 


0.61 ±0.22 


0.67±0.22 


0.39±0.22 


0.33±0.22 


1384 


0.0541 


7 


240±55 


0.58 


0.023 


0.71±0.26 


1.00±0.26 


0.29±0.26 


0.00±0.26 


3581 


0.0591 


5 


236±57 


0.57 


0.022 


0.50±0.25 


0.38±0.25 


0.79±0.20 


1.00±0.20 


3407 


0.0652 


4 


233±107 


0.56 


0.021 


0.29±0.26 


0.00±0.26 


0.29±0.26 


0.00±0.26 


2150 


0.0415 


9 


221 ±66 


0.53 


0.018 


0.50±0.34 


0.13±0.34 


0.50±0.34 


0.71 ±0.34 


2112 


0.0510 


12 


212±59 


0.51 


0.016 


0.71 ±0.26 


1.00±0.26 


0.29±0.26 


0.00±0.26 


3553 


0.0434 


10 


208±53 


0.50 


0.015 


0.61 ±0.22 


0.67±0.22 


0.61 ±0.22 


0.63±0.22 


3195 


0.0599 


11 


206±31 


0.50 


0.015 


0.69±0.19 


0.72±0.19 


0.50±0.20 


0.60±0.21 




U.U/40 


1 


ZU2± / / 


f\ A Q 


U.U14 


U.jU±U.2U 


0.47±0.20 


0.31±0.19 


f\ j_ A 1 O 

V.Z / ±\).ly 


1344 


0.0747 


1 


202±47 


0.48 


0.014 


0.77±0.15 


0.84±0.15 


0.50±0.18 


0.56±0.18 


3428 


0.0730 


5 


193+41 


0.46 


0.012 


0.39±0.22 


0.34±0.22 


0.61 ±0.22 


0.65±0.22 


3345 


0.0516 


4 


142±57 


0.34 


0.005 


0.71±0.26 


1.00±0.26 


0.29±0.26 


0.00±0.26 


3265 


0.0597 


4 


142±11 


0.34 


0.005 


0.50±0.25 


0.43±0.25 


0.50±0.25 


0.43±0.25 


3355 


0.0677 


4 


125+61 


0.30 


0.003 


0.84±0.16 


1.00±0.16 


0.16±0.16 


0.00±0.16 


3244 


0.0600 


4 


120±20 


0.29 


0.003 


0.71 ±0.26 


1.00±0.26 


0.71 ±0.26 


1.00±0.26 



3 From lvon der Linderi] i2007h 

4 Number of spectroscopically-confirmed cluster members from lvon der Lindenl d2007l) 



